Personalized pain management and anesthesia: preemptive risk identification and therapeutic decision support

ABSTRACT

Methods and compositions disclosed herein generally relate to methods of improving clinical and economic outcomes to address adverse effects related to anesthesia, analgesics, opioids, and inadequate pain relief. Embodiments of the invention relate to the association between genes, specific polymorphisms of genes, and non-genetic factors with inadequate pain relief and anesthesia-, analgesic, and/or opioid-related adverse effects. Embodiments of the invention can be used to determine and manage patient risk factors for development of adverse perioperative effects and can allow for personalized anesthesia and pain management for improvement of pain control and reduction of anesthesia-, analgesic-, and opioid-related adverse outcomes. These methods and compositions apply to non-surgical pain management with opioids. Therefore, patients who are genetically predisposed to risk of inadequate pain relief and/or serious side effects from anesthesia, analgesics, and/or opioids can be identified and individualized treatment plans developed for implementation by the clinician to improve clinical and economic outcomes.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application claims the benefit of priority under 35 U.S.C. § 119(e) to U.S. Provisional Application No. 61/565,400, METHODS AND COMPOSITIONS FOR CONDUCTING GENETIC PROFILING AND PERSONALIZATION OF SURGICAL ANALGESIA, filed on Nov. 30, 2011, and U.S. Provisional Application No. 61/661,073, PERSONALIZING PERIOPERATIVE OPIOID ANALGESIA-PREDICTIVE MODELS AND DECISION ALGORITHMS, filed on Jun. 18, 2012, which are incorporated by reference in their entirety.

STATEMENT REGARDING FEDERALLY-SPONSORED RESEARCH

This invention was made with U.S. Government support on behalf of the Department of Anesthesia, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio, USA (the Place Outcomes Research Award, Translational Research Initiative Grant, and Innovation Pilot Grant) and the National Center for Research Resources of the National Institutes of Health (NIH), USPHS Grant #UL1 RR026314, NIH Grant #eMERGE-II U01, and NIH eMERGE PGx supplement. The U.S. Government has certain rights in this invention.

FIELD OF THE INVENTION

The invention disclosed herein generally relates to methods and compositions comprising genetic markers for modulating surgical pain, anesthesia, post-surgical pain relief, and opioid-related adverse effects in a patient.

BACKGROUND

Safe and effective analgesia is an important medical and economic problem (Caldas, J. et al., Paediatr. Anaesth., 14:910-5 (2004); Duedahl, T. and Hansen, E., Paediatr. Anaesth., 17:756-74 (2007)). Approximately 28 million anesthetized surgical procedures are performed each year in the United States, many patients of which experience serious side-effects related to anesthesia, pain medications, and inadequate surgical pain relief. A significant fraction of the approximately 5 million children who undergo a painful surgery in the US each year experience inadequate pain relief and serious opioid-related side-effects (Cepeda, M., et al., Clin. Pharmacol. Ther., 74:102-12 (2003); Sadhasivam S., et al., Pediatrics, 129:832-8 (2012); Esclamado, R., et al., The Laryngoscope, 99:1125-9 (1989)). Safe and effective analgesia is an important unmet medical need, and its continued existence is an important clinical and perioperative safety and economic problem.

Adverse effects are observed throughout all classes of anesthetic treatments due to the narrow therapeutic indices of anesthetic and opioid pain medications. Moreover, a high degree of inter-individual variability in drug response underscores the challenges inherent to anesthetic treatment. Morphine, the most commonly used perioperative opioid, has a narrow therapeutic index and large inter-patient variations in analgesic response and serious side effects. Frequent inter-individual variations in responses to morphine have significant clinical and economic impact, with inadequate pain relief at one end of the spectrum of responses and serious adverse effects such as respiratory depression at the other end.

This inter-individual variability to drug response is presumed to be the result of a complex interaction of multiple factors. However, this complexity has stymied efforts to elucidate how genetic variability affects surgical pain and adverse responses to opioids. Accordingly, efforts to predict patient response to anesthesia can have a great impact in enabling clinicians to personalize analgesia to maximize pain relief while minimizing its adverse effects.

SUMMARY OF THE INVENTION

Embodiments of the invention encompass methods of determining a risk or a susceptibility to developing one or more adverse effects from administered anesthesia, analgesic, and/or opioid in a subject, the method including: obtaining a nucleic acid sample from a subject; analyzing the sample for presence or absence of at least one allele of at least one polymorphism associated with pain perception, persistent postoperative/chronic pain, and/or one or more anesthetic-, analgesic-, and/or opioid-related adverse effects; evaluating the subject for presence or absence of at least one non-genetic risk factor for suffering one or more adverse effects from administered anesthesia, analgesic, and/or opioid; determining the subject's risk or susceptibility to developing one or more adverse effects from administered anesthesia, analgesic, and/or opioid, wherein the presence of at least one non-genetic risk factor and/or the presence of absence of the at least one allele of the at least one polymorphism indicates that the patient has an elevated risk for suffering an one or more adverse anesthesia-, analgesic-, and/or opioid-related effects.

In some embodiments of the methods, the at least one non-genetic risk factor can be, for example, race, sex, age, weight, body mass index (BMI), or obstructive sleep apnea disease status. In some embodiments, the at least one non-genetic risk factor can be, for example, sleep disordered breathing, allergy status, sensitivity, a medical condition that affects an administrative route for delivered anesthetic, analgesic, and/or opioid, ethnicity, medical history, a drug interaction, psychological anxiety, stress level, or lifestyle.

In some embodiments of the methods, the at least one polymorphism associated with pain perception, persistent postoperative/chronic pain and/or one or more anesthetic adverse effects can be, for example, ABCB1 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-211, C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs 8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, V158M rs4680, the CYP2D6 *2, *2A, *3, *4,*5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19, *20, *35, *40, *41, *42, and *44 alleles, the PGRN-Seq genetic markers (Table 16), or markers in linkage disequilibrium therewith having r² value of at least 0.1 and/or |D′| value of at least 0.2.

In some embodiments, the sample can be analyzed for presence or absence of at least two polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or one or more anesthetic-, analgesic-, and/or opioid-related adverse effects. In some embodiments, the sample can be analyzed for the presence or absence of at least three polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or one or more anesthetic adverse effects. In some embodiments, the sample can be analyzed for the presence or absence of at least four polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or one or more anesthetic adverse effects.

In some embodiments, the subject can be a pediatric patient. In some embodiments, the subject can be at risk of inadequate pain relief and/or one or more serious side effects from administered anesthesia, analgesic, and/or opioids.

Embodiments of the invention also encompass compositions for identifying a patient at risk for suffering one or more adverse effects from administered anesthesia, analgesic, and/or opioid, including at least one polymorphism associated with pain perception, persistent postoperative/chronic pain, and/or one or more anesthetic-, analgesic-, and/or opioid-related adverse effects.

In some embodiments, the composition includes a gene chip. In some embodiments, the gene chip includes a low density array. In some embodiments, the composition includes a surface with a DNA array.

Embodiments of the invention also encompass methods of determining a risk or a susceptibility to developing obstructive sleep apnea (OSA) and/or a risk of OSA patients having increased pain and/or increased anesthetic-, analgesic-, and/or opioid-related adverse effects, the method including: obtaining a nucleic acid sample from a subject; analyzing the sample for the presence or absence of at least one allele of at least one polymorphism associated with pain perception and/or anesthetic-, analgesic, and/or opioid-related adverse effects; evaluating the subject for the presence or absence of at least one non-genetic risk factor for suffering adverse effects from administered anesthesia, analgesic, and/or opioid; determining the subject's risk of suffering adverse effects from administered anesthesia, analgesic, and/or opioid, wherein the presence of at least one non-genetic risk factor and/or the presence of absence of the at least one allele of the at least one polymorphism indicates that the patient has an elevated risk for susceptibility to developing OSA and/or a risk of OSA patients having increased pain and/or increased anesthetic-, analgesic-, and/or opioid-related adverse effects.

Embodiments of the invention also encompass methods of carrying out personalized treatment to maximize pain relief while minimizing adverse effects in a subject, the method including: identifying a subject in need of pain relief; obtaining a nucleic acid sample from the subject; analyzing the sample for the presence or absence of at least one allele of at least one polymorphism associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic-, analgesic, and/or opioid-related adverse effects; evaluating the subject for the presence or absence of at least one non-genetic risk factor for suffering adverse effects from administered anesthesia, analgesic, and/or opioid; determining the subject's risk of suffering adverse effects from administered anesthesia, analgesic, or opioid, wherein the presence of at least one non-genetic risk factor and/or the presence of absence of the at least one allele of the at least one polymorphism indicates that the patient has an elevated risk for suffering an adverse anesthesia-, analgesic-, and/or opioid-related effect; selecting one or more preferred therapeutics and/or alternatives and preferred dosage ranges based upon the determination; avoiding high-risk therapeutics based upon the determination; administering to the subject the one or more preferred therapeutics in the preferred dosage ranges, wherein administration of the one or more preferred therapeutics results in personalized treatment by maximizing pain relief while minimizing adverse effects.

In some embodiments, the therapeutic can be an analgesic, anesthetic, and/or opioid. In some embodiments, the analgesic, anesthetic, and/or opioid can be, for example, morphine or codeine. In some embodiments, the analgesic, anesthetic, and/or opioid can be, for example, benzocaine, chloroprocaine, cocaine, cyclomethycaine, dimethocaine/loracaine, piperocaine, propoxycaine, procaine/novocaine, proparacaine, tetracaine/amethocaine, articaine, bupivacaine, cinchocaine/dibucaine, etidocaine, levobupivacaine, lidocaine/lignocaine, mepivacaine, prilocaine, ropivacaine, trimecaine, saxitoxin, tetrodotoxin, menthol, eugenol, desflurane, enflurane, halothane, isoflurane, methoxyflurane, nitrous oxide, sevoflurane, xenon, amobarbital, methohexital, thiamylal, thiopental, diazepam, lorazepam, midazolam, etomidate, ketamine, propofol, alfentanil, fentanyl, remifentanil, sufentanil, buprenorphine, butorphanol, diamorphine, hydromorphone, levorphanol, meperidine, methadone, nalbuphine, oxycodone, hydrocodone, tramadol, ondansetron, dexamethasone, oxymorphone, pentazocine, succinylcholine, decamethonium, mivacurium, rapacuronium, atracurium, cisatracurium, rocuronium, vecuronium, alcuronium, doxacurium, gallamine, metocurine, pancuronium, pipecuronium, tubocurarine, morphine, fentanyl, hydromorphone, oxycodone, methadone, alfentanil, remifentanil, acetaminophen, NSAIDs, dexmedetomidine, clonidine, or combinations and derivatives thereof.

In some embodiments of the methods, the personalized treatment to maximize pain relief while minimizing adverse effects can be developed preemptively, prior to and/or during procedures and/or clinical care that involve administration of anesthetic, analgesic, and/or opioid. In some embodiments, the personalized treatment to maximize pain relief while minimizing adverse effects can be carried out postoperatively and in non-surgical pain management clinical settings.

BRIEF DESCRIPTION OF THE DRAWINGS

Those of skill in the art will understand that the drawings, described below, are for illustrative purposes only. The drawings are not intended to limit the scope of the present teachings in any way.

FIG. 1 depicts an exemplary hierarchical mixtures of experts (HME) decision tree, wherein the terminal nodes are called “experts” (or expert networks) that provide individual (and context dependent) prediction of the response.

FIG. 2 depicts a schematic showing the genetic and non-genetic factors linked to inter-individual variation in opioid adverse effects and analgesia.

FIG. 3 depicts a graph illustrating the odds ratio of ABCB1 genotypes and respiratory depression in a post-anesthesia care unit (PACU) stay.

FIG. 4 depicts graphs showing the odds ratio of ABCB1 and FAAH genotypes and post-operative nausea and vomiting (PONV) in a PACU stay.

FIGS. 5A-5D depict scatterplots and graphs illustrating the varying respiratory responses of ABCB1 genotypes to 5% carbon dioxide. FIG. 5A, before morphine administration; FIG. 5B, after morphine administration; FIG. 5C, resting minute ventilation, before and after surgery; and FIG. 5D, after CO₂ challenge, before and after surgery.

FIG. 6 depicts a plot showing the relationship between gene-gene (ABCB1 and FAAH) interactions and the probability of prolonged opioid-related respiratory depression.

FIGS. 7A-7D depict genetic variations of FAAH and hypercapnic ventilatory response in children and demonstrates that postoperative minute ventilation and ventilatory response to carbon dioxide were significantly lower compared to preoperative baseline values. Minute ventilation reduction differed between FAAH genotypes (p=0.007). FIG. 7A, before surgery and before CO₂ challenge; FIG. 7B after surgery but before CO₂ challenge; FIG. 7C before surgery but after CO₂ challenge; and FIG. 7D after surgery and after CO₂ challenge.

FIG. 8 depicts a plot illustrating the racial differences and predicted probabilities of postoperative morphine-induced adverse effects in children.

FIGS. 9A-9B depict the higher incidence of adverse effects experienced by Caucasian girls. FIG. 9A depicts the significant increase in PONV with increasing perioperative doses of morphine among Caucasian girls compared to Caucasian boys (p=0.0005); Caucasian boys were resistant to PONV and paradoxically had lower probability of PONV with increasing morphine doses. FIG. 9B depicts the significant increase in prolonged PACU stays due to opioid-related adverse effects among Caucasian girls compared to Caucasian boys (p=0.0113).

FIG. 10 depicts a scatterplot comparing morphine clearance in CCHMC data of children (6-15 years of age) and previously published data from Leiden (less than 3 years of age).

FIG. 11 depicts a scatterplot showing the clearance of morphine in African-American children (stars) and Caucasian children (diamonds).

FIG. 12 depicts a scatterplot illustrating the relationship between genetic markers and population stratification.

FIG. 13 depicts weight-normalized morphine clearance for OCT1 genotypes.

FIG. 14 depicts intravenous weight-normalized morphine clearance for OCT1 genotypes.

FIGS. 15A-15B depict the race difference in allele frequency and weight-normalized morphine clearance of OCT1 genotypes for African-Americans (15A) and Caucasians (15B).

FIG. 16 depicts morphine clearance for OCT genotypes in different weight groups.

FIG. 17 depicts morphine clearance with defective homozygous OCT1 genotypes.

FIG. 18 depicts the contribution of the rs622342 single nucleotide polymorphism (SNP) to the weight-normalized morphine clearance of OCT1 genotypes.

FIG. 19 depicts the weight and age dependency of the UGT2B7 polymorphism contribution.

FIG. 20 depicts acute opioid tolerance and the morphine dose requirement trends over 48 postoperative hours by TRPA genotype.

FIG. 21 depicts an “efficacy” decision tree for the post-operative analgesic (PA) use vs. no intervention (noi) classification, using SNP data and covariates such as race. Leaves (nodes) of the tree that represent strata with increased risk of inadequate pain relief and those with relatively lower risk are indicated on the figure; the number of patients in each of the two classes in each node is shown explicitly.

FIG. 22 depicts a “safety” decision tree for the classification of adverse effects (AE) vs. the other classes (NOS). Leaves corresponding to strata with increased risk of AE are indicated on the figure. Only the two main leaves that can be classified most easily by combinations of ABCB1, ADRB2, COMT, and FAAH polymorphisms are shown explicitly within the middle branch.

FIG. 23 depicts a decision tree for the ability to predict obstructive sleep apnea (OSA) vs the other classes (NS), using the established panel of SNPs. The OSA tree has about 72% prediction accuracy (with about 50% baseline).

FIG. 24 depicts patterns of linkage disequilibrium (LD) and haplotype blocks around the FAAH gene in Caucasian (CEU) and African ancestry in the southwest US (ASW) HapMap 3 cohorts. Standard D prime measure is used, with dark squares corresponding to high LD (for reference, see <http colon slash slash>www<dot>hapmap<dot>org).

FIG. 25 depicts the patient cohort for the study of codeine adverse drug reactions (ADRs) and analgesic effects.

FIG. 26 depicts the CYP2D6 total activity score (TAS) association with the predicted number of ADRs.

FIGS. 27A-27B depict the factors contributing to the probability of sedation after codeine administration. FIG. 27A, pain score=30; FIG. 27B, pain score=60. POD, post-operative days.

FIGS. 28A-28B depict the factors contributing to pain intensity before and after codeine administration in 7 year olds (28A) and in 12 year olds (28B).

DETAILED DESCRIPTION OF THE INVENTION Definitions

Unless otherwise noted, terms are to be understood according to conventional usage by those of ordinary skill in the relevant art.

A “single nucleotide polymorphism,” or “SNP,” is a DNA sequence variation occurring when a single nucleotide at a specific location in the genome differs between members of a species or between paired chromosomes in an individual. Most SNP polymorphisms have two alleles. Each individual is in this instance either homozygous for one allele of the polymorphism (i.e. both chromosomal copies of the individual have the same nucleotide at the SNP location), or the individual is heterozygous (i.e. the two sister chromosomes of the individual contain different nucleotides). The SNP nomenclature as reported herein refers to the official Reference SNP (rs) ID identification tag as assigned to each unique SNP by the National Center for Biotechnological Information (NCBI) or identifies the residue change associated with the identified polymorphism. SNP genotyping arrays have become an important tool for cohort identification and stratification, phenotype-genotype association studies, discovery of disease markers, prediction of molecular phenotypes, and clinical decision support.

As used herein, the term “variant” refers to a segment of DNA that differs from the reference DNA.

As used herein, the term “marker” or a “polymorphic marker” refers to a variant. Alleles that differ from the reference are referred to as “variant” alleles.

A “microsatellite” is a polymorphic marker that has multiple small repeats of bases that are typically 2-8 nucleotides in length (such as CA repeats) at a particular site, in which the number of repeat lengths varies in the general population. An “indel” is a common form of polymorphism comprising a small insertion or deletion that is typically only a few nucleotides long.

As used herein, the term “haplotype” refers to a segment of genomic DNA that is characterized by a specific combination of a series of polymorphic markers arranged along the segment. For diploid organisms such as humans, a haplotype comprises one member of the pair of alleles for each polymorphic marker or locus along the segment. In some embodiments, the haplotype can comprise an allele for each of two or more markers, three or more markers, four or more markers, or five or more markers. As used herein, the term “susceptibility” refers to the proneness of an individual towards the development of a certain state (e.g. a certain trait, phenotype, or disease), or towards being less able to resist a particular state than the average individual. The term encompasses both increased susceptibility and decreased susceptibility. Thus, particular alleles at polymorphic markers and/or haplotypes of the invention as described herein can be characteristic of increased susceptibility (i.e. increased risk) of adverse anesthetic effects, as characterized by a relative risk (RR) or odds ratio (OR) of greater than one for the particular allele or haplotype. Alternatively, some markers and/or haplotypes of the invention can be characteristic of decreased susceptibility (i.e. decreased risk) of adverse anesthetic effects, as characterized by a relative risk of less than one.

As used herein, the term “and/or” shall in the present context be understood to indicate that either or both of the items connected by it are involved. In other words, the term herein shall be taken to mean “one or the other or both.”

As used herein, the term “look-up table” is a table that correlates one form of data to another form, or one or more forms of data to a predicted outcome to which the data is relevant, such as a phenotype or trait. For example, a look-up table can comprise a correlation between allelic data for at least one polymorphic marker and a particular trait or phenotype, such as a particular disease diagnosis, that an individual who comprises the particular allelic data is likely to display, or is more likely to display than individuals who do not comprise the particular allelic data. Look-up tables can be multidimensional, i.e. they can contain information about multiple alleles for single markers simultaneously, or they can contain information about multiple markers; they can also comprise other factors, such as particulars about diseases diagnoses, racial information, biomarkers, biochemical measurements, therapeutic methods, or drugs, etc.

A “nucleic acid sample” as described herein, refers to a sample obtained from an individual that contains nucleic acid (DNA or RNA). In certain embodiments, such as, for example, the detection of specific polymorphic markers and/or haplotypes, the nucleic acid sample comprises genomic DNA. Such a nucleic acid sample can be obtained from any source that contains genomic DNA, including a blood sample, sample of amniotic fluid, sample of cerebrospinal fluid, or tissue sample from skin, muscle, buccal or conjunctival mucosa, placenta, gastrointestinal tract, or other organs.

As used herein, the term “decision tree” refers to a standard machine learning technique for multivariate data analysis and classification (Hastie, T., et al. The Elements of Statistical Learning, Second Edition, Springer (2009); Witten, I. and Frank, E. Data Mining: Practical Machine Learning Tools and Techniques, Second Edition, Morgan Kaufmann, San Francisco (2005); Hothorn, T. Journal of Computational and Graphical Statistics, 15:651-74 (2010)). Decision trees can be used to derive easily interpretable and intuitive rules for decision support systems.

“Therapeutic” and “therapeutics” as used herein refers primarily to one or more anesthetic, analgesic, and/or opioid compositions having an actual or potential beneficial effect for a patient. In some embodiments, the term can also include accompanying modes of treatment and/or administration and/or co-administration with other compositions and/or treatments, as recognized by those of skill in the art of anesthesia and analgesia.

“Alternative” and “alternatives” as used herein refers to one or more compositions and/or treatment and/or administration and/or co-administration that is generally considered by those of skill in the art to be an alternative to a given therapeutic.

A “high risk therapeutic” as used herein refers to a therapeutic and/or alternative that, in comparison to other therapeutics and/or alternatives, implicates an elevated risk of adverse effects, side effects, and/or other negative risks or eventualities.

Inter-individual variability to drug response is a significant clinical and perioperative safety and economic problem. Embodiments of the invention described herein include methods and compositions for carrying out personalized analgesia treatment to maximize pain relief while minimizing its adverse effects for patients.

Approximately 50% of the inter-individual variability in clinical response to morphine can be explained by SNPs in the genes involved in pain mechanisms and opioid metabolism, transport, and receptor signaling. Though previous research has identified individual genes contributing to the pain processing pathway by looking at single genotype-phenotype associations, no work has heretofore investigated how the interplay of multiple genes, gene-gene interactions, and non-genetic factors, including demographic, phenotypic and other data, affect analgesic response, nor has any research demonstrated the ability to predict a patient's level of pain upon treatment, or to facilitate both basic and translational clinical research using such information.

Embodiments of the invention are directed to the identification and management of risk factors of adverse postoperative outcomes, and personalized anesthesia, opioid administration, and pain management for improvement of pain control, analgesia, and reduction of anesthesia- and opioid-related adverse outcomes. This can include examining the influence of specific polymorphisms of genes involved in pain perception, opioid transport, and liver drug metabolism, amongst others, on postoperative pain and the analgesia and adverse effects of opioids and anesthesia. Steps and benchmarks can include: 1. Identifying patients who are genetically predisposed to risk of inadequate pain relief and/or serious side effects from anesthesia and opioids including morphine. 2. Exploring the effects of race and associated genetic profile on postoperative pain control and adverse effects with morphine and other opioids in patients. 3. Developing and validating a predictive and personalized clinical decision model to improve clinical, perioperative, and economic outcomes.

Embodiments of the invention are also directed to methods of identifying a patient at risk for suffering adverse effects from administered anesthesia, analgesic, or opioid, comprising obtaining a sample from the patient and analyzing the sample for the presence or absence of at least one polymorphism associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic adverse effects, wherein the presence of absence of the at least one polymorphism indicates that the patient has an elevated risk for suffering an adverse opioid- and/or anesthesia-related effect.

Embodiments of the invention are also directed to methods of identifying a patient at risk for suffering from obstructive sleep apnea (OSA), comprising obtaining a sample from the patient and analyzing the sample for the presence or absence of at least one polymorphism associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic adverse effects and/or OSA, wherein the presence of absence of the at least one polymorphism indicates that the patient has an elevated risk for suffering from OSA.

In some embodiments, the sample is analyzed for the presence or absence of at least two polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic-, analgesic-, and/or opioid-related adverse effects. In some embodiments, the sample is analyzed for the presence or absence of at least three polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic-, analgesic-, and/or opioid-related adverse effects. In some embodiments, the sample is analyzed for the presence or absence of at least four polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic-, analgesic-, and/or opioid-related adverse effects.

As disclosed herein, particular marker alleles or haplotypes (e.g. ABCB1 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-211C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, V158M rs4680, the CYP2D6 *2, *2A, *3, *4,*5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19, *20, *35, *40, *41, *42, and *44 alleles, the PharmacoGenomics Research Network (PGRN) genetic markers on the PGRN-Seq genetic markers platform (Table 16), and markers in linkage disequilibrium therewith can be associated with development of adverse effects to one or more administered anesthesia, analgesic, and/or opioid. In some embodiments, the marker allele or haplotype is one that confers a significant risk or susceptibility to developing such effects. Accordingly, in embodiments of the invention, a method of determining a risk or a susceptibility to developing adverse effects from administered anesthesia, analgesic, and/or opioid in a subject is provided, the method comprising determining the presence or absence of at least one allele of at least one polymorphic marker in a nucleic acid sample obtained from the subject, wherein the at least one polymorphic marker is selected from the group consisting of: ABCB1 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-211C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, V158M rs4680, the CYP2D6 *2, *2A, *3, *4,*5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19, *20, *35, *40, *41, *42, and *44 alleles, the PGRN-Seq genetic markers (Table 16), and markers in linkage disequilibrium therewith.

In some embodiments, a method of determining a risk or a susceptibility to developing adverse effects from administered anesthesia, analgesic, and/or opioid in a subject is provided, the method comprising screening for at least one marker selected from ABCB1 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-211C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, the CYP2D6 *2, *2A, *3, *4,*5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19, *20, *35, *40, *41, *42, and *44 alleles, the PGRN-Seq genetic markers (Table 16), and markers in linkage disequilibrium therewith. In some embodiments, the marker allele or haplotype is more frequently present in a subject having, or who is susceptible to, developing adverse effects from administered anesthesia (affected), as compared to the frequency of its presence in a healthy subject (control, such as population controls). In some embodiments, the significance of association of the at least one marker allele or haplotype is characterized by a p-value <0.05. In some embodiments, the significance of association is characterized by smaller p-values, such as <0.01, <0.001, <0.0001, <0.00001, <0.000001, <0.0000001, <0.00000001 or <0.000000001.

Embodiments of the invention are also directed to a composition for identifying a patient at risk for suffering adverse effects from administered anesthesia, analgesic, and/or opioid, comprising at least one polymorphism associated with pain perception, persistent postoperative/chronic pain, and/or analgesic and anesthetic adverse effects. In some embodiments, the composition comprises at least two polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic adverse effects. In some embodiments, the composition comprises at least three polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic adverse effects. In some embodiments, the composition comprises at least four polymorphisms associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic adverse effects. In some embodiments, the composition comprises a gene chip. In some embodiments, the composition comprises a surface with a DNA array.

Embodiments of the invention include use of each gene or each polymorphism separately as a diagnostic or prognostic marker, or use of a few marker genes combined into a panel display format so that several marker genes can be detected to increase reliability and efficiency. Further, any of the genes or polymorphisms identified as disclosed herein can be used individually or as a set of genes in any combination with any of the other genes or other polymorphisms that are disclosed in the application.

Exemplary Markers

In relatively small adult studies, specific polymorphisms of the catechol-o-methyltransferase (COMT), mu opioid receptor (OPRM1), and the ATP binding cassette B1 (ABCB1) (also known as MCR1) genes have been associated with high pain sensitivity, increased morphine requirement and side effects. As disclosed herein, heretofore unknown genetic variations were identified that represent high risk factors of opioid-related adverse effects and poor pain control. These genetic variations and risk factors are predictive of potential outcomes and can be used to create decision trees and algorithms based upon these biomarkers to enable personalized interventions that can be clinically implemented. These decision trees and algorithms can be used to evaluate the efficacy and safety of administering opioids and/or analgesics.

The use of the decision trees and algorithms described herein can improve perioperative clinical outcomes by providing better analgesia with minimal adverse effects, as well as economic outcomes, by reducing prolonged PACU stays, emergency room visits for inadequate pain control, dehydration following opioid-induced vomiting, and respiratory depression requiring oxygen and intense respiratory support measures. These clinical decision algorithms can enable tailored opioid selection and dosing to maximize pain relief while minimizing serious adverse effects. These findings can be extrapolated to adults and children with surgical or non-surgical pain requiring opioids.

The results described herein disclose critically important findings regarding inter-individual variations in pain perception and opioid responses. This information provides the essential foundation for clinical decision algorithms that enable tailored opioid selection and dosing so as to maximize pain relief while minimizing serious adverse-effects in children. Millions of children have the potential to benefit from this research, which will improve safety, quality of care, and economics of perioperative care.

Drug transporters at the blood-brain barrier, including the ATP binding cassette B1 (ABCB1) transporter protein, have the potential to significantly influence the clinical efficacy and safety of opioids. SNPs can therefore impact inter-individual variability in response to opioids, and gene variants have been shown to modulate cerebral pharmacokinetics of morphine and increase its analgesic and side effects. For example, adult carriers of the GG variant of the mu opioid receptor (OPRM1) SNP, A118G, can require up to a 2-4 fold higher dose of morphine than AA variants.

As disclosed herein, preliminary experiments comparing genotypes of ABCB1 (rs1045642) from 150 children undergoing tonsillectomy showed 4 to 10 fold inter-individual variations in morphine's respiratory depression and post-operative nausea and vomiting (PONV). Furthermore, allelic frequency of the TT genotype in ABCB1, which correlates with a higher risk of respiratory depression and PONV, was 27%.

Fatty acid amide hydrolase (FAAH) is part of the endocannabinoid system, a neuromodulator of vomiting. FAAH degrades anandamide, as an endogenous agonist for CB1 receptors. Canniboid agonists are effective antiemetics against opioids through CB1 receptors. As disclosed herein, genetic variations in FAAH (rs4141964) can have strong associations with opioid adverse effects, especially PONV and prolonged PACU stays. As further disclosed herein, genetic variations in TRPA1 (rs1443952) can also have strong associations with opioid adverse effects.

Catechol-o-methyltransferase (COMT) is a key regulator of pain perception, cognitive function, and affective mood. COMT is also a pivotal regulator of catecholamines concentrations in the pain perception pathway. Haplotypes of four common COMT SNPs (rs6269, rs4633, rs4818, and rs4680) have been shown to be associated with three levels of pain perceptions (low, average, and high pain sensitivity). Additionally, the V158M polymorphism (rs4680) influences the human response to pain. As disclosed herein, there can be strong associations between COMT and inter-individual variations in pain scores and opioid requirements. Certain GCHI and other genotypes are associated with persistent postoperative/chronic pain.

OSA is over-represented among African-American patients, who in turn are more likely to be at a higher risk of inadequate pain relief, as indicated by the efficacy tree described herein. OSA and race therefore represent important covariates; this is further highlighted by the finding that ADRB2 and FAAH gene polymorphisms can be used to relatively accurately predict race (with ˜80% prediction accuracy). However, OSA itself seems to have a relatively strong, although poorly understood, genetic component (Larkin, E. et al. Resp. Crit. Care Med 182:947-53 (2010)). As disclosed herein, OSA can be predicted using SNPs in just two genes (ADRB2 rs1042717 and ABCB1). A simple decision rule was developed that expresses the risk of OSA; this decision rule achieves about 73% classification accuracy (compared to 50% for a baseline classifier).

As described herein, a systematic study was conducted to identify much less understood genetic variants underlying clinical responses to opioids in children with the goal of improving and personalizing the postoperative care and pain management in children. A standardized site study protocol allowed for unambiguous assignment and accurate quantification of the observed phenotypes. Genotype data was collected using a specialized SNP panel. The candidate genes included ABCB1, COMT, OPRM1, FAAH, ADRB2 and a number of other genes that were chosen based on their allele frequencies and clinical evidences of important associations in adults with opioid analgesic and adverse effects.

Preliminary statistical analyses revealed several significant associations between genetic, as well as non-genetic factors and postoperative opioid adverse effects and inadequate analgesia. In particular, the TT genotype of the ABCB1 SNP rs1045642 (C3435T) was found (after adjusting for obstructive sleep apnea) to be associated with a higher risk of morphine-induced respiratory depression than the CC genotype. In the ABCB1 TT genotype, resting minute ventilation (MV) after morphine decreased by 47.5% compared to only 18.4% in CC and CT genotypes (p<0.05). At the same time, a number of other, relatively weak associations with SNPs in FAAH, COMT, and other genes, as well as indication of epistatic interactions between, e.g., ABCB1 and FAAH (as well as ABCB1 and ADRB2) SNPs were found.

Different strata with distinct patterns of such interactions (with ABCB1 playing a prominent role in some strata) were identified. These data also suggest that African American children had inadequate pain control, while Caucasian children had higher incidence of adverse effects from similar doses of morphine.

Significant differences were found in the allelic frequency of ABCB1 (and other genes). For example, the TT genotype of ABCB1 SNP rs1045642 that predisposes children to opioid induced respiratory depression (with >4 fold higher incidence than the CC genotype) was found to be as frequent as 27% in Caucasian children, as compared to 2-3% in our African-American children.

In some embodiments, the at least one polymorphism that is analyzed is one that is found in a gene from selected from the group of: ABCB1, ABCC2, ABCC3, ABCDI, ADRB2, ANKK1, beta-arrestin-2, COMT, CHRM3, CNRI, CNR2, CYP2D6, CYP3A4, CYP2C9, DRD2, eCNOS, FAAH, FLJ37396, GCHI, GRIN2B, HTR2A, ILRN, IL2RB, KCNB2, KCNJ6, LOC100287246, MDR1, MC1R, NBPF22P, NTRK1, OCT1, OPRM1, PTGSI, PTGS2, SLC6A2, SLC6A4, SCN9A, Stat6, TANK, TRPA1, TRPV1, TNF-alpha, TACR1, UGT2B7, UGT2B14, UGT2B15, UGT2B17, UGT1A1, UGT1A3, UGT1A6, UGT1A8, V158M, and ZNF429.

In some embodiments, the at least one polymorphism that is analyzed is selected from the group of: ABCB1 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-21 IC>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, V158M rs4680 activity, the CYP2D6 *2, *2A, *3, *4, *5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19, *20, *35, *40, *41, *42, and *44 alleles, and the PGRN-Seq genetic markers (Table 16).

The polymorphisms as disclosed herein can be correlated with pain perception, persistent postoperative/chronic pain, anesthetic or opioid requirement, and anesthetic or opioid adverse effects. In some embodiments, the polymorphisms as disclosed herein are correlated with pain perception, persistent postoperative/chronic pain, anesthetic or opioid requirement, and anesthetic-, analgesic-, and/or opioid-related adverse effects in a pediatric patient. In some embodiments, the polymorphisms as disclosed herein can be correlated with OSA.

Assessment for Markers and Haplotypes

The genomic sequence within populations is not identical when individuals are compared. Rather, the genome exhibits sequence variability between individuals at many locations in the genome. Such variations in sequence are commonly referred to as polymorphisms, and there are many such sites within each genome. For example, the human genome exhibits sequence variations which occur on average every 500 base pairs. The most common sequence variant consists of base variations at a single base position in the genome, and such sequence variants, or polymorphisms, are commonly called SNPs. These SNPs are believed to have occurred in a single mutational event, and therefore there are usually two possible alleles possible at each SNP site: the original allele and the mutated allele. Due to natural genetic drift and possibly also selective pressure, the original mutation has resulted in a polymorphism characterized by a particular frequency of its alleles in any given population. Many other types of sequence variants are found in the human genome, including mini- and microsatellites, and insertions, deletions, and inversions (also called copy number variations (CNVs)). A polymorphic microsatellite has multiple small repeats of bases (such as CA repeats, TG on the complimentary strand) at a particular site in which the number of repeat lengths varies in the general population. In general terms, each version of the sequence with respect to the polymorphic site represents a specific allele of the polymorphic site. These sequence variants can all be referred to as polymorphisms, occurring at specific polymorphic sites characteristic of the sequence variant in question. In general terms, polymorphisms can comprise any number of specific alleles. Thus, in some embodiments of the invention, the polymorphism is characterized by the presence of two or more alleles in any given population. In some embodiments, the polymorphism is characterized by the presence of three or more alleles. In some embodiments, the polymorphism is characterized by four or more alleles, five or more alleles, six or more alleles, seven or more alleles, nine or more alleles, or ten or more alleles. All such polymorphisms can be utilized in the methods and kits as disclosed herein.

Reference can be made to different alleles at a polymorphic site without choosing a reference allele. Alternatively, a reference sequence can be referred to for a particular polymorphic site. The reference allele can be referred to as the “wild-type” allele, and it usually is chosen as either the first sequenced allele or as the allele from a “non-affected” individual (e.g. an individual that does not display a trait or disease phenotype).

Alleles for SNP markers as referred to herein refer to the bases A, C, G, or T as they occur at the polymorphic site in the SNP assay employed. The allele codes for SNPs used herein are as follows: I=A, 2-C, 3=G, 4=T. However, the person skilled in the art will understand that by assaying or reading the opposite DNA strand, the complementary allele can in each case be measured. Thus, for a polymorphic site (polymorphic marker) characterized by an A/G polymorphism, the assay employed can be designed to specifically detect the presence of one or both of the two bases possible, i.e. A and G. Alternatively, by designing an assay that is designed to detect the complimentary strand on the DNA template, the presence of the complementary bases T and C can be measured. Quantitatively (for example, in terms of relative risk), identical results would be obtained from measurement of either DNA strand (+ strand or − strand).

Typically, a reference sequence is referred to for a particular sequence. Alleles that differ from the reference are sometimes referred to as “variant” alleles. A variant sequence, as used herein, refers to a sequence that differs from the reference sequence but is otherwise substantially similar. Alleles at the polymorphic genetic markers described herein are variants. Variants can include changes that affect a polypeptide. Sequence differences, when compared to a reference nucleotide sequence, can include the insertion or deletion of a single nucleotide, or of more than one nucleotide, resulting in a frame shift; the change of at least one nucleotide, resulting in a change in the encoded amino acid; the change of at least one nucleotide, resulting in the generation of a premature stop codon; the deletion of several nucleotides, resulting in a deletion of one or more amino acids encoded by the nucleotides; the insertion of one or several nucleotides, such as by unequal recombination or gene conversion, resulting in an interruption of the coding sequence of a reading frame; duplication of all or a part of a sequence; transposition; or a rearrangement of a nucleotide sequence. Such sequence changes can alter the polypeptide encoded by the nucleic acid. For example, if the change in the nucleic acid sequence causes a frame shift, the frame shift can result in a change in the encoded amino acids, and/or can result in the generation of a premature stop codon, causing generation of a truncated polypeptide. Alternatively, a polymorphism associated with a disease or trait can be a synonymous change in one or more nucleotides (i.e. a change that does not result in a change in the amino acid sequence). Such a polymorphism can, for example, alter splice sites, affect the stability or transport of mRNA, or otherwise affect the transcription or translation of an encoded polypeptide. It can also alter DNA to increase the possibility that structural changes, such as amplifications or deletions, occur at the somatic level. The polypeptide encoded by the reference nucleotide sequence is the “reference” polypeptide with a particular reference amino acid sequence, and polypeptides encoded by variant alleles are referred to as “variant” polypeptides with variant amino acid sequences.

A haplotype refers to a segment of DNA that is characterized by a specific combination of alleles arranged along the segment. For diploid organisms such as humans, a haplotype comprises one member of the pair of alleles for each polymorphic marker or locus. In a certain embodiment, the haplotype can comprise two or more alleles, three or more alleles, four or more alleles, or five or more alleles, each allele corresponding to a specific polymorphic marker along the segment. Haplotypes can comprise a combination of various polymorphic markers, e.g. SNPs and microsatellites, having particular alleles at the polymorphic sites. The haplotypes thus comprise a combination of alleles at various genetic markers.

Detecting specific polymorphic markers and/or haplotypes can be accomplished by methods known in the art for detecting sequences at polymorphic sites. For example, standard techniques for genotyping for the presence of SNPs and/or microsatellite markers can be used, such as fluorescence-based techniques (e.g. Chen, X. et al., Genome Res. 9(5): 492-98 (1999); Kutyavin et al., Nucleic Acid Res. 34:e128 (2006)), utilizing PCR, LCR, Nested PCR, and other techniques for nucleic acid amplification. Specific commercial methodologies available for SNP genotyping include, but are not limited to, TaqMan genotyping assays and SNPlex platforms (Applied Biosystems), gel electrophoresis (Applied Biosystems), mass spectrometry (e.g. MassARRAY system from Sequenom), minisequencing methods, real-time PCR, Bio-Plex system (BioRad), CEQ and SNPstream systems (Beckman), array hybridization technology (e.g. Affymetrix GeneChip; Perlegen), BeadArray Technologies (e.g. Illumina GoldenGate and Infinium assays), array tag technology (e.g. Parallele), and endonuclease-based fluorescence hybridization technology (Invader, Third Wave). Thus, by use of these or other methods available to the person skilled in the art, one or more alleles at polymorphic markers, including microsatellites, SNPs, or other types of polymorphic markers, can be identified.

As disclosed herein, an individual who is at an increased susceptibility (i.e. increased risk) for a disease or trait is an individual in whom at least one specific allele at one or more polymorphic marker or haplotype conferring increased susceptibility (increased risk) for the disease or trait is identified (i.e. at-risk marker alleles or haplotypes). The at-risk marker or haplotype is one that confers an increased risk (increased susceptibility) of the disease. In some embodiments, significance associated with a marker or haplotype is measured by a relative risk (RR). In some embodiments, significance associated with a marker or haplotype is measured by an odds ratio (OR). In a further embodiment, the significance is measured by a percentage. In some embodiments, a significant increased risk is measured as a risk (relative risk and/or odds ratio) of at least 1.2, including but not limited to: at least 1.2, at least 1.3, at least 1.4, at least 1.5, at least 1.6, at least 1.7, 1.8, at least 1.9, at least 2.0, at least 2.5, at least 3.0, at least 4.0, and at least 5.0. In some embodiments, a risk (relative risk and/or odds ratio) of at least 1.2 is significant. In some embodiments, a risk of at least 1.3 is significant. In some embodiments, a risk of at least 1.4 is significant. In some embodiments, a relative risk of at least 1.5 is significant. In some embodiments, a significant increase in risk is at least 1.7 is significant. However, other cutoffs are also contemplated, e.g. at least 1.15, 1.25, 1.35, and so on, and such cutoffs are also within scope of the present invention. In some embodiments, a significant increase in risk is at least about 20%, including but not limited to about 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 100%, 150%, 200%, 300%, and 500%. In some embodiments, a significant increase in risk is at least 20%. In some embodiments, a significant increase in risk is at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, and at least 100%. Other cutoffs or ranges as deemed suitable by the person skilled in the art to characterize the invention are however also contemplated, and those are also within scope of the present invention. In some embodiments, a significant increase in risk is characterized by a p-value, such as a p-value of less than 0.05, less than 0.01, less than 0.001, less than 0.0001, less than 0.00001, less than 0.000001, less than 0.0000001, less than 0.00000001, or less than 0.000000001.

An at-risk polymorphic marker or haplotype as described herein is one where at least one allele of at least one marker or haplotype is more frequently present in an individual at risk for the disease (or trait) (affected), or diagnosed with the disease, compared to the frequency of its presence in a comparison group (control), such that the presence of the marker or haplotype is indicative of susceptibility to the disease. In some embodiments, the control group can be a population sample, i.e. a random sample from the general population. In some embodiments, the control group is represented by a group of individuals who are disease-free. In some embodiments, such disease-free controls can be characterized by the absence of one or more specific disease-associated symptoms. Alternatively, the disease-free controls are those that have not been diagnosed with the disease. In some embodiments, the disease-free control group is characterized by the absence of one or more disease-specific risk factors. In some embodiments, such risk factors include at least one environmental risk factor. Representative environmental factors are natural products, minerals, or other chemicals which are known to affect, or contemplated to affect, the risk of developing the specific disease or trait. Other environmental risk factors are risk factors related to lifestyle, including but not limited to food and drink habits, geographical location of main habitat, and occupational risk factors. In some embodiments, the risk factors comprise at least one additional genetic risk factor.

An example of a simple test for correlation would be a Fisher-exact test on a two by two table. Given a cohort of chromosomes, the two by two table is constructed out of the number of chromosomes that include both of the markers or haplotypes, one of the markers or haplotypes but not the other, and neither of the markers or haplotypes. Other statistical tests of association known to the skilled person are also contemplated and are also within scope of the invention.

In some embodiments of the invention, an individual who is at a decreased susceptibility (i.e. at a decreased risk) for a disease or trait is an individual in whom at least one specific allele at one or more polymorphic marker or haplotype conferring decreased susceptibility for the disease or trait is identified. The marker alleles and/or haplotypes conferring decreased risk are also said to be protective. In some embodiments, the protective marker or haplotype is one that confers a significant decreased risk (or susceptibility) of the disease or trait. In some embodiments, significant decreased risk is measured as a relative risk (or odds ratio) of less than 0.9, including but not limited to less than 0.9, less than 0.8, less than 0.7, less than 0.6, less than 0.5, less than 0.4, less than 0.3, less than 0.2, and less than 0.1. In some embodiments, significant decreased risk is less than 0.7. In some embodiments, significant decreased risk is less than 0.5. In some embodiments, significant decreased risk is less than 0.3. In some embodiments, the decrease in risk (or susceptibility) is at least 20%, including but not limited to at least 25%, at least 30%, at least 35%, at least 40%, at least 45%, at least 50%, at least 55%, at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 90%, at least 95%, and at least 98%. In some embodiments, a significant decrease in risk is at least about 30%. In some embodiments, a significant decrease in risk is at least about 50%. In some embodiments, the decrease in risk is at least about 70%. Other cutoffs or ranges as deemed suitable by the person skilled in the art to characterize the invention are however also contemplated, and those are also within scope of embodiments of the invention.

The person skilled in the art will appreciate that for markers with two alleles present in the population being studied (such as SNPs), and wherein one allele is found in increased frequency in a group of individuals with a trait or disease in the population, compared with controls, the other allele of the marker will be found in decreased frequency in the group of individuals with the trait or disease, compared with controls. In such a case, one allele of the marker (the one found in increased frequency in individuals with the trait or disease) will be the at-risk allele, while the other allele will be a protective allele.

A genetic variant associated with a disease or a trait can be used alone to predict the risk of the disease for a given genotype. For a biallelic marker, such as an SNP, there are three possible genotypes: homozygote for the at-risk variant, heterozygote, and non-carrier of the at-risk variant. Risk associated with variants at multiple loci can be used to estimate overall risk. For multiple SNP variants, there are k possible genotypes, with k=3^(n)×2^(P); where n is the number autosomal loci, and p the number of gonosomal (sex chromosomal) loci. Overall risk assessment calculations for a plurality of risk variants typically assume that the relative risks of different genetic variants multiply, i.e. the overall risk (e.g. RR or OR) associated with a particular genotype combination is the product of the risk values for the genotype at each locus. If the risk presented is the relative risk for a person or a specific genotype for a person, compared to a reference population with matched gender and ethnicity, then the combined risk is the product of the locus-specific risk values and which also corresponds to an overall risk estimate compared with the population. If the risk for a person is based on a comparison to non-carriers of the at-risk allele, then the combined risk corresponds to an estimate that compares the person with a given combination of genotypes at all loci to a group of individuals who do not carry risk variants at any of those loci. The group of non-carriers of any at risk variant has the lowest estimated risk and has a combined risk, compared with itself (i.e. non-carriers) of 1.0 but has an overall risk, compared with the population, of less than 1.0. It should be noted that the group of non-carriers can potentially be very small, especially for a large number of loci, and, in that case, its relevance is correspondingly small.

The multiplicative model is a parsimonious model that usually fits the data of complex traits reasonably well. Deviations from multiplicity have been rarely described in the context of common variants for common diseases and if reported are usually only suggestive, since very large sample sizes are typically required to be able to demonstrate statistical interactions between loci.

By way of example, consider a total of eight variants that have been described to associate with prostate cancer (Gudmundsson, et al., Nat Genet 39:631-7 (2007); Gudmundsson, et al., Nat Genet 39:977-83 (2007); Yeager, et al., Nat Genet 39:645-49 (2007); Amundadottir, et al., Nat Genet 38:652-8 (2006); Haiman, et al., Nat Genet 39:638-44 (2007)). Seven of these loci are on autosomes, and the remaining locus is on chromosome X. The total number of theoretical genotypic combinations is then 3⁷×2¹=4374. Some of those genotypic classes are very rare but are still possible, and these can be considered for overall risk assessment. The multiplicative model applied in the case of multiple genetic variants can also be valid in conjugation with non-genetic risk variants, assuming that the genetic variant does not clearly correlate with the “environmental” factor. In other words, genetic and non-genetic at-risk variants can be assessed under the multiplicative model to estimate combined risk, assuming that the non-genetic and genetic risk factors do not interact.

Using the same quantitative approach, the combined or overall risk associated with a plurality of variants associated with adverse anesthetic effects can be assessed, including combinations of any one of the markers and/or polymorphisms as disclosed herein or markers in linkage disequilibrium therewith.

Linkage Disequilibrium

The natural phenomenon of recombination, which occurs on average once for each chromosomal pair during each meiotic event, represents one way in which nature provides variations in sequence (and biological function by consequence). It has been discovered that recombination does not occur randomly in the genome; rather, there are large variations in the frequency of recombination rates, resulting in small regions of high recombination frequency (also called recombination hotspots) and larger regions of low recombination frequency, which are commonly referred to as linkage disequilibrium (LD) blocks (Myers, et al., Biochem Soc Trans 34:526-30 (2006); Jeffreys, et al., Nat Genet 29:217-22 (2001); May, et al., Nat Genet 31:272-5 (2002)).

Linkage disequilibrium (LD) refers to a non-random assortment of two genetic elements. For example, if a particular genetic element (e.g. an allele of a polymorphic marker, or a haplotype) occurs in a population at a frequency of 0.50 (50%), and another element occurs at a frequency of 0.50 (50%), then the predicted occurrence of a person's having both elements is 0.25 (25%), assuming a random distribution of the elements. However, if it is discovered that the two elements occur together at a frequency higher than 0.25, then the elements are said to be in linkage disequilibrium, since they tend to be inherited together at a higher rate than what their independent frequencies of occurrence (e.g. allele or haplotype frequencies) would predict. Roughly speaking, LD is generally correlated with the frequency of recombination events between the two elements. Allele or haplotype frequencies can be determined in a population by genotyping individuals in a population and determining the frequency of the occurrence of each allele or haplotype in the population. For populations of diploids, e.g. human populations, individuals will typically have two alleles or allelic combinations for each genetic element (e.g. a marker, haplotype, or gene).

Many different measures have been proposed for assessing the strength of linkage disequilibrium (reviewed in Devlin, B. & Risch, N., Genomics 29:311-22 (1995)). Most capture the strength of association between pairs of bi-allelic sites. Two important pairwise measures of LD are r² (sometimes denoted A²) and |D′| (Lewontin, R., Genetics 49:49-67 (1964); Hill, W. G. & Robertson, A. Theor. Appl. Genet. 22:226-231 (1968)). Both measures range from 0 (no disequilibrium) to 1 (“complete” disequilibrium), but their interpretation is slightly different. |D′| is defined in such a way that it is equal to 1 if just two or three of the possible haplotypes are present, and it is <1 if all four possible haplotypes are present. Therefore, a value of |D′| that is <1 indicates that historical recombination can have occurred between two sites (recurrent mutation can also cause |D′| to be <1, but for SNPs, this is typically regarded as being less likely than recombination). The measure r² represents the statistical correlation between two sites and takes the value of 1 if only two haplotypes are present.

The r² measure is a relevant measure for association mapping because there is a simple inverse relationship between r² and the sample size sufficient to detect association between susceptibility loci and SNPs. These measures are defined for pairs of sites, but, for some applications, a determination of how strong LD is across an entire region that contains many polymorphic sites can be desirable (e.g. testing whether the strength of LD differs significantly among loci or across populations, or whether there is more or less LD in a region than predicted under a particular model). Measuring LD across a region is not straightforward, but one approach is to use the measure r, which was developed in population genetics. Roughly speaking, r measures how much recombination would be sufficient under a particular population model to generate the LD that is seen in the data. This type of method can also provide a statistically rigorous approach to the problem of determining whether LD data provide evidence for the presence of recombination hotspots. For the methods described herein, a significant r² value can be at least 0.1, such as at least 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.91, 0.92, 0.93, 0.94, 0.95, 0.96, 0.97, 0.98, or at least 0.99. In some embodiments, the significant r² value can be at least 0.2. Alternatively, linkage disequilibrium as described herein refers to linkage disequilibrium characterized by |D′| values of at least 0.2, such as 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.85, 0.9, 0.95, 0.96, 0.97, 0.98, or at least 0.99. Thus, linkage disequilibrium represents a correlation between alleles of distinct markers. It is measured by correlation coefficient or |D′| (r² up to 1.0 and |D′| up to 1.0). In some embodiments, linkage disequilibrium is defined in terms of values for both the r² and |D′| measures. In some embodiments, a significant linkage disequilibrium is defined as r²>0.1 and |D′|>0.8. In some embodiments, a significant linkage disequilibrium is defined as r²>0.2 and |D′|>0.9. Other combinations and permutations of values of r² and |D′| for determining linkage disequilibrium are also contemplated, and are also within the scope of embodiments of the invention. Linkage disequilibrium can be determined in a single human population, as defined herein, or it can be determined in a collection of samples comprising individuals from more than one human population. In some embodiments, LD is determined in a sample from one or more of the HapMap populations (Caucasian, African, Japanese, Chinese) (as defined at http<colon slash slash>www<dot>hapmap<dot>org).

If all polymorphisms in the genome were independent at the population level (i.e. no LD), then every single one of them would need to be investigated in association studies, to assess all the different polymorphic states. However, due to linkage disequilibrium between polymorphisms, tightly linked polymorphisms are strongly correlated, which reduces the number of polymorphisms that need to be investigated in an association study to observe a significant association. Another consequence of LD is that many polymorphisms can give an association signal due to the fact that these polymorphisms are strongly correlated.

Genomic LD maps have been generated across the genome, and such LD maps have been proposed to serve as framework for mapping disease-genes (Risch, N. & Merkiangas, K, Science 273:1516-1517 (1996); Maniatis, et al., Proc Natl Acad Sci USA 99:2228-2233 (2002); Reich, et al., Nature 411:199-204 (2001)).

It is now established that many portions of the human genome can be broken into series of discrete haplotype blocks containing a few common haplotypes; for these blocks, linkage disequilibrium data provides little evidence indicating recombination (see, e.g., Wall., J. D. and Pritchard, J. K., Nature Reviews Genetics 4:587-597 (2003); Daly, et al., Nature Genet. 29:229-232 (2001); Gabriel, et al., Science 296:2225-2229 (2002); Patil, et al., Science 294:1719-1723 (2001); Dawson, et al., Nature 418:544-548 (2002); Phillips, et al., Nature Genet. 33:382-387 (2003)).

There are two main methods for defining these haplotype blocks: blocks can be defined as regions of DNA that have limited haplotype diversity (see, e.g., Daly, et al., Nature Genet. 29:229-232 (2001); Patil, et al., Science 294:1719-1723 (2001); Dawson, et al., Nature 418:544-548 (2002); Zhang, et al., Proc. Natl. Acad Sci. USA 99:7335-7339 (2002)), or as regions between transition zones having extensive historical recombination, identified using linkage disequilibrium (see, e.g., Gabriel, et al., Science 296:2225-2229 (2002); Phillips, et al., Nature Genet. 33:382-387 (2003); Wang, et al., Am. J. Hum. Genet. 71:1227-1234 (2002); Stumpf, M. P., and Goldstein, D. B., Curr. Biol. 13:1-8 (2003)). More recently, a fine-scale map of recombination rates and corresponding hotspots across the human genome has been generated (Myers, et al., Science 310:321-32324 (2005); Myers, et al., Biochem Soc Trans 34:526530 (2006)). The map reveals the enormous variation in recombination across the genome, with recombination rates as high as 10-60 cM/Mb in hotspots, while closer to 0 in intervening regions, which thus represent regions of limited haplotype diversity and high LD. The map can therefore be used to define haplotype blocks/LD blocks as regions flanked by recombination hotspots. As used herein, the terms “haplotype block” or “LD block” includes blocks defined by any of the above described characteristics, or other alternative methods used by the person skilled in the art to define such regions.

Haplotype blocks (LD blocks) can be used to map associations between phenotype and haplotype status, using single markers or haplotypes comprising a plurality of markers. The main haplotypes can be identified in each haplotype block, and a set of “tagging” SNPs or markers (the smallest set of SNPs or markers sufficient to distinguish among the haplotypes) can then be identified. These tagging SNPs or markers can then be used in assessment of samples from groups of individuals, in order to identify association between the phenotype and haplotype. If desired, neighboring haplotype blocks can be assessed concurrently, as there can also exist linkage disequilibrium among the haplotype blocks.

It has thus become apparent that for any given observed association to a polymorphic marker in the genome, it is likely that additional markers in the genome also show association. This is a natural consequence of the uneven distribution of LD across the genome, as observed by the large variation in recombination rates. The markers used to detect association thus in a sense represent “tags” for a genomic region (i.e., a haplotype block or LD block) that is associating with a given disease or trait and as such are useful for use in the methods and kits of the present invention. One or more causative (functional) variants or mutations can reside within the region found to be associating to the disease or trait. The functional variant can be another SNP, a tandem repeat polymorphism (such as a minisatellite or a microsatellite), a transposable element, or a copy number variation, such as an inversion, deletion, or insertion. Such variants in LD with the variants described herein can confer a higher relative risk (RR) or odds ratio (OR) than observed for the tagging markers used to detect the association. The present invention thus refers to the markers used for detecting association to the disease, as described herein, as well as markers in linkage disequilibrium with the markers. Thus, in some embodiments of the invention, markers that are in LD with the markers and/or haplotypes of the invention, as described herein, can be used as surrogate markers. The surrogate markers have in some embodiments relative risk (RR) and/or odds ratio (OR) values smaller than for the markers or haplotypes initially found to be associating with the disease, as described herein. In some embodiments, the surrogate markers have RR or OR values greater than those initially determined for the markers initially found to be associating with the disease, as described herein. An example of such an embodiment would be a rare or relatively rare (such as <10% allelic population frequency) variant in LD with a more common variant (>10% population frequency) initially found to be associating with the disease, such as the variants described herein. Identifying and using such markers for detecting the association discovered by the inventors as described herein can be performed by routine methods well-known to the person skilled in the art and are therefore within the scope of the present invention.

Determination of Haplotype Frequency

The frequencies of haplotypes in patient and control groups can be estimated using an expectation-maximization algorithm (Dempster, et al., J. R. Stat. Soc. 8, 39:1-38 (1977)). An implementation of this algorithm that can handle missing genotypes and uncertainty with the phase can be used. Under the null hypothesis, the patients and the controls are assumed to have identical frequencies. Using a likelihood approach, an alternative hypothesis is tested, where a candidate at-risk haplotype, which can include the markers described herein, is allowed to have a higher frequency in patients than controls, while the ratios of the frequencies of other haplotypes are assumed to be the same in both groups. Likelihoods are maximized separately under both hypotheses, and a corresponding 1-df likelihood ratio statistic is used to evaluate the statistical significance.

To look for at-risk and protective markers and haplotypes within a susceptibility region, for example within an LD block, association of all possible combinations of genotyped markers within the region is studied. The combined patient and control groups can be randomly divided into two sets, equal in size to the original group of patients and controls. The marker and haplotype analysis is then repeated, and the most significant p-value registered is determined. This randomization scheme can be repeated, for example, over 100 times to construct an empirical distribution of p-values. In some embodiments, a p-value of <0.05 is indicative of a significant marker and/or haplotype association.

Risk Assessment and Diagnostics

Within any given population, there is an absolute risk of developing a disease or trait, defined as the chance of a person developing the specific disease or trait over a specified time period. For example, a woman's lifetime absolute risk of breast cancer is one in nine. That is to say, one woman in every nine will develop breast cancer at some point in their lives. Risk is typically measured by looking at very large numbers of people, rather than at a particular individual. Risk is often presented in terms of absolute risk (AR) and relative risk (RR). Relative risk is used to compare risks associated with two variants or the risks of two different groups of people. For example, it can be used to compare a group of people with a certain genotype with another group having a different genotype. For a disease or trait, a relative risk of 2 means that one group has twice the chance of developing the disease or trait as the other group. The risk presented is usually the relative risk for a person, or a specific genotype of a person, compared to the population with matched gender and ethnicity. Risks of two individuals of the same gender and ethnicity can be compared in a simple manner. For example, if, compared to the population, the first individual has a relative risk of 1.5, and the second has a relative risk of 0.5, then the risk of the first individual compared to the second individual is 1.5/0.5=3.

Risk Calculations

The creation of a model to calculate the overall genetic risk involves two steps: i) conversion of odds-ratios for a single genetic variant into relative risk, and ii) combination of risk from multiple variants in different genetic loci into a single relative risk value.

Deriving Risk from Odds-Ratios

Most gene discovery studies for complex diseases or traits that have been published to date in authoritative journals have employed a case-control design due to their retrospective setup. These studies sample and genotype a selected set of cases (people who have the specified disease condition) and control individuals. The interest is in genetic variants (alleles) whose frequency in cases and controls differ significantly.

The results are typically reported in odds ratios, which describe the ratio between the fraction (probability) with the risk variant (carriers) versus the non-risk variant (non-carriers) in the groups of affected versus the controls, i.e. expressed in terms of probabilities conditional on the affection status:

OR=(Pr(c|A)/Pr(nc|A))/(Pr(c|C)/Pr(nc|C))

In some embodiments, the absolute risk for the disease or trait is what is determined, i.e. the fraction of those individuals carrying the risk variant who get the disease or, in other words, the probability of getting the disease. This number cannot be directly measured in case-control studies, in part because the ratio of cases versus controls is typically not the same as that in the general population. However, under certain assumptions, the risk can be calculated from the odds ratio value.

Combining the Risk from Multiple Markers

When genotypes of many SNP variants are used to estimate the risk for an individual, unless otherwise stated, a multiplicative model for risk can be assumed. This means that the combined genetic risk relative to the population is calculated as the product of the corresponding estimates for individual markers, e.g. for two markers g1 and g2:

RR(g1,g2)=RR(g1)RR(g2)

The underlying assumption is that the risk factors occur and behave independently, i.e. that the joint conditional probabilities can be represented as products:

Pr(A|g1,g2)=Pr(A|g1)Pr(A|g2)/Pr(A) and Pr(g1,g2)=Pr(g1)Pr(g2)

In embodiments where markers are closely spaced on the genome, i.e. are in linkage disequilibrium such that the concurrence of two or more risk alleles is correlated, a haplotype modeling where the odds-ratios are defined for all allele combinations of the correlated SNPs can be employed.

As an example, consider an individual who has the following genotypes at four markers associated with risk of type-2 diabetes along with the risk relative to the population at each marker.

Chromo 3 PPARG CC Calculated risk: RR(CC)=1.03

Chromo 6 CDKAL1 GG Calculated risk: RR(GG)=1.30

Chromo 9 CDKN2A AG Calculated risk: RR(AG)=0.88

Chromo 11 TCF7L2 TI Calculated risk: RR(TI)=1.54

Combined, the overall risk relative to the population for this individual is: 1.03×1.30×0.88×1.54=1.81.

Risk Assessment for Adverse Anesthetic Effects

As disclosed herein, certain polymorphic markers and haplotypes comprising such markers are found to be useful for risk assessment of adverse reactions to surgical anesthesia. Risk assessment can involve the use of the markers for determining a susceptibility to such adverse effects. Particular alleles of polymorphic markers (e.g. SNPs) are found more frequently in individuals with particular susceptibility to such adverse effects. Therefore, these marker alleles have predictive value for determining whether these individuals will suffer post-operative adverse effects from administered anesthesia, or a susceptibility to the same. Tagging markers in linkage disequilibrium with at-risk variants (or protective variants) as disclosed herein can be used as surrogates for these markers (and/or haplotypes). Such surrogate markers can be located within a particular haplotype block or LD block. Such surrogate markers can also sometimes be located outside the physical boundaries of such a haplotype block or LD block, either in close vicinity of the LD block/haplotype block or possibly also located in a more distant genomic location.

Long-distance LD can arise if particular genomic regions (e.g. genes) are in a functional relationship. For example, if two genes encode proteins that play a role in a shared metabolic pathway, then particular variants in one gene can have a direct impact on observed variants for the other gene. An example is a case where a variant in one gene leads to increased expression of the gene product. To counteract this effect and preserve overall flux of the particular pathway, this variant can have led to selection of one (or more) variants at a second gene that confers decreased expression levels of that gene. These two genes can be located in different genomic locations, possibly on different chromosomes, but variants within the genes are in apparent LD, not because of their shared physical location within a region of high LD but rather due to evolutionary forces. Such LD is also contemplated and within the scope of embodiments of the invention. The skilled person will appreciate that many other scenarios of functional gene-gene interaction are possible, and the particular example discussed here represents only one such possible scenario.

Markers with values of r² equal to 1 are equivalent surrogates for the at-risk variants, i.e. genotypes for one marker perfectly predict genotypes for the other. Markers with smaller values of r² than 1 can also be surrogates for the at-risk variant or alternatively represent variants with relative risk values as high as or possibly even higher than the at-risk variant. The at-risk variant identified is not necessarily the functional variant itself but is in this instance in linkage disequilibrium with the true functional variant. The functional variant can, for example, be a tandem repeat, such as, but not limited to, a minisatellite or a microsatellite, a transposable element (e.g. an A/u element), or a structural alteration, such as a deletion, insertion, or inversion (sometimes also called copy number variations, or CNVs). Embodiments of the invention encompass the assessment of such surrogate markers for the markers as disclosed herein. Such markers are annotated, mapped, and listed in public databases, as is well-known to the skilled person, or can alternatively be readily identified by sequencing the region or a part of the region identified by the markers of the present invention in a group of individuals and identify polymorphisms in the resulting group of sequences. As a consequence, the person skilled in the art can readily and without undue experimentation genotype surrogate markers in linkage disequilibrium with the markers and/or haplotypes as described herein. The tagging or surrogate markers in LD with the detected at-risk variants also have predictive value for detecting association to adverse anesthetic effects, or a susceptibility to the same, in an individual. These tagging or surrogate markers that are in LD with the identified markers as disclosed herein can also include other markers that distinguish among haplotypes, as these similarly have predictive value for detecting susceptibility to the particular trait.

Embodiments of the invention can be directed to methods of assessing a sample comprising genomic DNA from an individual for the presence of variants described herein to be associated with adverse anesthetic effects. Such methods typically comprise steps that detect the presence or absence of at least one allele of at least one polymorphic marker, using methods well-known to the skilled person and further described herein, and, based on the outcome of such assessment, determine whether the individual from whom the sample is derived is at increased or decreased risk (increased or decreased susceptibility) of adverse anesthetic effects. In some embodiments, detecting particular alleles of polymorphic markers can be carried out by obtaining nucleic acid sequence data about a particular human individual that identifies at least one allele of at least one polymorphic marker. Different alleles of the at least one marker are associated with different susceptibility to such effects in humans. Obtaining nucleic acid sequence data can comprise identifying the nucleic acid sequence at a single nucleotide position, which is sufficient to identify alleles at SNPs. The nucleic acid sequence data can also comprise sequence information at any other number of nucleotide positions, in particular for genetic markers that comprise multiple nucleotide positions, and can be anywhere from two to hundreds of thousands, possibly even millions, of nucleotides (in particular, in the case of CNVs).

In some embodiments, the methods comprise utilization of a dataset comprising information about the genotype status of at least one polymorphic marker associated with a disease or trait (or markers in linkage disequilibrium with at least one marker associated with the disease or trait). In other words, a dataset containing information about such genetic status, for example in the form of sequence data, genotype counts at a certain polymorphic marker, or a plurality of markers (e.g. an indication of the presence or absence of certain at-risk alleles), or actual genotypes for one or more markers, can be queried for the presence or absence of certain at-risk alleles at certain polymorphic markers as disclosed herein to be associated with the disease or trait. A positive result for a variant (e.g. marker allele) associated with the disease or trait is indicative of the individual from whom the dataset is derived is at increased susceptibility (increased risk) of the disease.

In some embodiments, a polymorphic marker is correlated to a disease or trait by referencing genotype data for the polymorphic marker to a look-up table that comprises correlations between at least one allele of the polymorphism and the disease. In some embodiments, the table comprises a correlation for one polymorphism. In some embodiments, the table comprises a correlation for a plurality of polymorphisms. By referencing to a look-up table that gives an indication of a correlation between a marker and the disease, a risk for the disease or trait or a susceptibility to the disease or trait can be identified in the individual from whom the sample is derived. In some embodiments, the correlation is reported as a statistical measure. The statistical measure can be reported as a risk measure, such as a relative risk (RR), an absolute risk (AR), or an odds ratio (OR).

The markers described herein can be useful for risk assessment and diagnostic purposes, either alone or in combination. The results of the risk of suffering adverse anesthetic effects based on the markers described herein can also be combined with data for other genetic markers or risk factors for adverse anesthetic effects to establish an overall risk. Thus, even in cases where the increase in risk by individual markers is relatively modest, e.g. on the order of 10-30%, the association can have significant implications. Thus, relatively common variants can have significant contribution to the overall risk (population attributable risk is high), or combinations of markers can be used to define groups of individuals who, based on the combined risk of the markers, are at significant combined risk of developing complications due to adverse anesthetic effects.

Thus, in some embodiments, a plurality of variants (genetic markers, biomarkers, and/or haplotypes) is used for overall risk assessment. In some embodiments, these variants can be selected from the variants as disclosed herein. In some embodiments, the use of the variants as disclosed herein in combination with other variants known to be useful for diagnosing a susceptibility to adverse anesthetic effects is used. In such embodiments, the genotype status of a plurality of markers and/or haplotypes is determined in an individual, and the status of the individual is compared with the population frequency of the associated variants or the frequency of the variants in clinically healthy subjects, such as age-matched and sex-matched subjects. Methods known in the art, such as multivariate analyses or joint risk analyses or other methods known to the skilled person, can subsequently be used to determine the overall risk conferred based on the genotype status at the multiple loci. Assessment of risk based on such analysis can subsequently be used in the methods, uses, and kits of the invention, as described herein.

Individuals who are homozygous for at-risk variants for suffering adverse anesthetic effects can be at particularly high risk of developing such effects. This is due to the dose-dependent effect of at-risk alleles, such that the risk for homozygous carriers is generally estimated as the risk for each allelic copy squared.

As disclosed herein, the haplotype block structure of the human genome has the effect that a large number of variants (markers and/or haplotypes) in linkage disequilibrium with the variant originally associated with a disease or trait can be used as surrogate markers for assessing association to the disease or trait. The number of such surrogate markers will depend on factors such as the historical recombination rate in the region, the mutational frequency in the region (i.e. the number of polymorphic sites or markers in the region), and the extent of LD (size of the LD block) in the region. These markers are usually located within the physical boundaries of the LD block or haplotype block in question as defined using the methods described herein, or by other methods known to the person skilled in the art. However, sometimes marker and haplotype association is found to extend beyond the physical boundaries of the haplotype block as defined, as disclosed herein. Such markers and/or haplotypes can in those cases be also used as surrogate markers and/or haplotypes for the markers and/or haplotypes physically residing within the haplotype block as discussed above. As a consequence, markers and haplotypes in LD (typically characterized by inter-marker r² values of greater than 0.1, such as r² greater than 0.2, including r² greater than 0.3, and also including markers correlated by values for r² greater than 0.4) with the markers and haplotypes of the present invention are also within the scope of the invention, even if they are physically located beyond the boundaries of the haplotype block as defined. This includes markers that are disclosed herein but can also include other markers that are in strong LD (e.g. characterized by r² greater than 0.1 or 0.2 and/or |D′|>0.8) with the described markers.

For the SNP markers described herein, the opposite allele to the allele found to be in excess in patients (at-risk allele) is found in decreased frequency in situations where adverse effects from anesthetic and analgesic agents are observed. These markers and haplotypes in LD and/or comprising such markers along with other factors, including age, race, weight (or body mass index, BMI), are thus protective for adverse anesthetic and analgesic effects, i.e. they confer a decreased risk or susceptibility of individuals carrying these markers and/or haplotypes developing pain and/or adverse effects from anesthetics and analgesics.

In some embodiments, variants as disclosed herein, including certain haplotypes, can comprise a combination of various genetic markers, e.g. SNPs and microsatellites. Detecting haplotypes can be accomplished by methods known in the art and/or described herein for detecting sequences at polymorphic sites. Furthermore, correlation between certain haplotypes or sets of markers and disease phenotypes can be verified using standard techniques. A representative example of a simple test for correlation would be a Fisher-exact test on a two by two table.

In some embodiments, a marker allele or haplotype found to be associated with adverse anesthetic effects is one in which the marker allele or haplotype is more frequently present in an individual at risk for suffering such effects (affected), compared to the frequency of its presence in a healthy individual (control), or in a randomly selected individual from the population, wherein the presence of the marker allele or haplotype is indicative of a susceptibility to suffering adverse anesthetic effects. In some embodiments, at-risk markers in linkage disequilibrium with one or more markers shown herein to be associated with suffering adverse anesthetic effects (e.g. marker alleles as listed in Tables 1 and 2 and Example 4) are tagging markers that are more frequently present in an individual at risk for suffering adverse anesthetic effects (affected), compared to the frequency of their presence in a healthy individual (control) or in a randomly selected individual from the population, wherein the presence of the tagging markers is indicative of increased susceptibility to suffering such effects. In some embodiments, at-risk marker alleles (i.e. conferring increased susceptibility) in linkage disequilibrium with one or more markers found to be associated with adverse anesthetic effects are markers comprising one or more allele that is more frequently present in an individual at risk for suffering such effects, compared to the frequency of their presence in a healthy individual (control), wherein the presence of the markers is indicative of increased susceptibility to suffering from adverse anesthetic effects.

Predictive Models

A decision tree is a standard machine learning technique for multivariate data analysis and classification that can be used to derive easily interpretable and intuitive rules for decision support systems. Decision tress can be viewed as a recursive partitioning approach, in which data is hierarchically divided into strata by simple logical rules. The advantage of decision trees is their simplicity, ability to handle categorical and numerical variables, as well as missing values, robustness to outliers and scaling, and the ability to combine feature selection with stratification and classification. Decision trees can also be used to derive easy to interpret and intuitive rules for decision support systems.

As described herein, decision trees are used to select and combine the most predictive SNPs with demographic, clinical, and other input features into simple logical rules that can be used to classify patients and predict adverse effects, thereby enabling robust and accurate point-of-care prediction of inadequate pain relief and opioid-related adverse effects. Such knowledge allows for individualized treatment.

As described herein, the observed distinct strata and complex interaction patterns have resulted in the systematic determination of whether interactions of specific polymorphisms of genes, such as those involved in opioid transport, sensing, and metabolism, significantly influence morphine's adverse effects in children. Due to the limitations of current approaches, there is a need for tailored solutions and extensions that can provide more robust and accurate decision rules for personalized interventions.

Distinct strata associated with specific patterns of gene-gene interactions in the context of adverse effects, pain sensitivity, and other clinical phenotypes were identified (with ABCB1 playing a prominent role in some strata) and analyzed first using standard recursive partitioning (or decision tree-based) approaches, such as CART. Genotyping, demographic, and relevant clinical data were then combined to derive logical rules for the prediction of patients with high vs. low risk of adverse effects, including opioid-induced respiratory depression outcomes. The predictive power of such models is limited by inherent biological noise, limited sample sizes, and complex pattern of interactions representing multiple mechanisms that can lead to variation in molecular and eventually clinical outcomes. Weaker associations that can manifest themselves only in the context of specific strata can be present as well but are more difficult to detect Decision trees and associated logical rules were enhanced by incorporating allelic (additive), dominant, and recessive models, as well as ternary trees wherein each genotype value decoupled were also implemented to facilitate identification and analysis of distinct strata. In addition, haplotype reconstruction for each gene (strong patterns of LD are observed for several genes included here) was incorporated, using Phase and similar population-based models, allowing for the simplification of decision rules and accounting for many implicit interactions observed in the data. Classical decision trees, support vector machines (SVMs), and hierarchical mixtures of experts (HME) models were integrated.

The HME approach can be viewed as a probabilistic decision tree (see Hasti et al., chapter 9 and references therein). An example of an HME decision tree is shown in FIG. 1, where the terminal nodes are called “experts” (or expert networks) that provide individual (and context-dependent) prediction of the response (in this case, adverse effects). These individual predictions are combined by the decision models (“gating networks”) of the non-terminal nodes. Typically, a linear or logistic regression model is fit in each terminal node, using a linear combination of inputs, as for example in the following model for objective postoperative pain score (P) as a function of genes and other predictors:

E(P|X)=β₀+β₁Age+β₂Sex+β₃Race+β₄GCH1+β₅COMT+β₆ABCB1+Other Factors

In the above equation, E(Y|X) denotes the expected value of Y given a vector of predictor variables (features) X. Expert models can be generalized using other suitable classification (or regression) approaches, such as by using robust linear SVM predictors. In this approach, terminal SVM models can be regarded as an ensemble of expert classifiers that are optimized in a context-dependent manner (including one-class SVMs when applicable).

This can be compared with developing an ensemble of classifiers optimized on different subsets of the data (such as in the context of cross-validation), which are combined to provide more robust results and associated confidence levels (which are defined in terms of consistency within the ensemble). Here, such subsets were defined by robust distinct strata obtained using standard decision trees (possibly with re-sampling and boosting as well). Thus, the tree topology of such an extended HME model would be expected to represent stable strata identified in the original decision tree analysis, addressing one of the limitations of the HME approach, namely the lack of methods for finding a good tree topology. Therefore, the resulting tree preserves the ease of interpretation of the decision rules, while improving accuracy.

Other standard statistical and machine learning methods, including neural networks, prototype, and kernel-based approaches, were also applied to further dissect patterns of gene-gene interactions and to assess their predictive power (and to evaluate relative merits of the new approach). Cross-validation, in which the overall training cohort is repeatedly (and randomly) divided into training and validation subsets, was used to evaluate the accuracy and to assess the stability of the resulting decision rules and observed strata. Standard measures of accuracy, including overall classification accuracy, precision, recall, and area under ROC curve, were used to evaluate the accuracy of the decision rules based on the proposed mixed recursive partitioning models.

As described herein, a systematic multivariate analysis of associations was initiated between gene-gene interactions and other confounding factors (such as race, age, weight or BMI, etc.) and the respective outcomes, using standard CART and C4.5 decision trees, as well as other machine learning approaches. These methods were first used to identify and analyze potential patterns of gene-gene interactions and other factors predictive of inadequate pain relief or pain sensitivity (represented by the need for the post-operative analgesic (PA) use).

The data described herein demonstrate that African-American children had inadequate pain control and Caucasian children had a higher incidence of adverse effects from similar doses of morphine (Sadhasivam S., et al., Pediatrics, 129:832-8(2012)). Concordant differences in allelic frequency of ABCB1 (and other genes) were observed. For example, the T genotype of ABCB1 SNP rs1045642 that predisposes children to opioid-induced respiratory depression (>4 fold higher incidence than CC genotype) was found with 27% frequency in Caucasian children, as compared to 2-3% in African-American children.

For example, the following exemplary rules were developed and applied to the decision trees described in Examples 23 and 26, respectively:

-   -   IF(Race=African-American) THEN high risk of pain sensitivity and         inadequate pain relief;     -   ELSE IF((Race=White) AND (GCHI rs441417=TIT)) THEN moderate risk         of pain sensitivity and inadequate pain relief;     -   ELSE IF((Race=White) AND (GCHI rs441417=CC OR CT)) THEN low risk         of pain sensitivity and inadequate pain relief.     -   and     -   IF(((ADRB2 rs1042717=GG) AND (MDR1 rs1045642-CC)) OR ((ADRB2         rs1042717=AA or GA) AND (MDR1 rs9282564=AA))) THEN increased         risk of OSA

Exemplary embodiments of the invention involve use of the above rules to stratify populations and provide a risk assessment for pain sensitivity and OSA risk.

Study Population

In a general sense, the methods and kits of the invention can be utilized from samples containing nucleic acid material (DNA or RNA) from any source and from any individual or from genotype data derived from such samples. In some embodiments, the individual is a human individual. The individual can be an adult, child, or fetus. The nucleic acid source can be any sample comprising nucleic acid material, including biological samples, or a sample comprising nucleic acid material derived therefrom. Embodiments of the invention also provide for assessing markers and/or haplotypes in individuals who are members of a target population. Such a target population is in one embodiment a population or group of individuals at risk of developing adverse anesthetic effects, based on other parameters such as, for example, genetic factors, biomarkers, biophysical parameters, history of anesthesia-related complications, allergic reactions to anesthesia, family history of anesthesia-related complications, and the like.

In some embodiments, the target population includes individuals from specific age subgroups, such as those under the age of 18, under the age of 16, or under the age of 12 or 10. The individuals can be of either sex, males or females.

It is believed that the markers found to be associated with adverse anesthetic, analgesia, and/or opioid effects as disclosed herein can show similar association in other human populations outside of the population employed in the current study. In some embodiments, the human subjects are from one or more human populations or ethnic groups, including, but not limited to, Caucasian populations, European populations, American populations, Eurasian populations, Asian populations, Central/South Asian populations, East Asian populations, Middle Eastern populations, African populations, Hispanic populations, and Oceanian populations.

The racial contribution in individual subjects can also be determined by genetic analysis. Genetic analysis of ancestry can be carried out using unlinked microsatellite markers such as those set out in Smith et al. (Am J Hum Genet 74, 1001-13 (2004)).

In some embodiments, the invention relates to markers and/or haplotypes identified in specific populations, as described in the above. The person skilled in the art will appreciate that measures of LD can give different results when applied to different populations. This is due to different population histories of different human populations as well as differential selective pressures that can have led to differences in LD in specific genomic regions. It is also well known to the person skilled in the art that certain markers, e.g. SNP markers, have different population frequencies in different populations or are polymorphic in one population but not in another. The person skilled in the art will however apply the methods available and as taught herein to practice embodiments of the invention in any given human population. This can include assessment of polymorphic markers in the LD region of the present invention, so as to identify those markers that give the strongest association within the specific population. Thus, the at-risk variants of the present invention can reside on different haplotype backgrounds and in different frequencies in various human populations. However, utilizing methods known in the art and the markers of the present invention, the invention can be practiced in any given human population.

All references cited herein are incorporated by reference in their entirety as though fully set forth. Unless defined otherwise, technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Singleton et al., Dictionary of Microbiology and Molecular Biology 3^(rd) ed., J. Wiley & Sons (New York, N.Y. 2001); March, Advanced Organic Chemistry Reactions, Mechanisms and Structure 5^(th) ed., J. Wiley & Sons (New York, N.Y. 2001); and Sambrook and Russel, Molecular Cloning: A Laboratory Manual 3^(rd) ed, Cold Spring Harbor Laboratory Press (Cold Spring Harbor, N.Y. 2001) provide one skilled in the art with a general guide to many of the terms used in the present application.

One skilled in the art will recognize many methods and materials similar or equivalent to those described herein, which can be used in the practice of the present invention. Indeed, the present invention is in no way limited to the methods and materials described.

Diagnostic Methods

In embodiments of the invention, methods of evaluating, or aiding in the evaluation of, a patient at risk for suffering adverse effects from administered anesthetic, analgesia, and/or opioid, are provided, the methods comprising detecting particular alleles at genetic markers that are correlated with such effects. In some embodiments, methods to determine susceptibility to developing adverse effects to administered anesthetic, analgesia, and/or opioid are provided, the methods comprising detecting at least one allele of at least one polymorphic marker (e.g. the markers described herein). As disclosed herein, particular alleles of particular markers or haplotypes are indicative of a susceptibility to suffering adverse effects from administered analgesia. Prognostic or predictive assays for evaluating a patient's risk or susceptibility to such effects can also be used to develop a customized anesthetic protocol for a subject prior to and during procedures that involve administration of anesthetic, analgesic, and/or opioid.

In some embodiments, methods of identifying a patient at risk for suffering adverse effects from administered anesthetic, analgesia, and/or opioid are carried out by analyzing samples for the presence of absence of at least one polymorphism associated with pain perception, persistent postoperative/chronic pain, and/or anesthetic-, analgesic-, and/or opioid-related adverse effects using microarray or gene chip technology, wherein the microarray or gene chip comprises the at least one polymorphism.

In some embodiments, the gene chip comprises a low density array.

In some embodiments, the methods include evaluation of additional clinical information to tailor pain therapy and anesthesia and minimize anesthetic and analgesic medication-related adverse effects. In some embodiments, the methods include evaluation of additional clinical information to tailor pain therapy and opioids and minimize opioid and analgesic medication-related adverse effects. Additional clinical information for use in such methods include, but are not limited to, patient age, race, weight (or BMI), sleep disordered breathing, gender, allergies, sensitivities, or medical conditions that affect administrative routes for delivered anesthetic, analgesia, and/or opioid. Other factors, such as, for example, ethnicity, medical history, drug interactions, psychological anxiety, stress level, and lifestyle can also be evaluated as part of the methods.

Embodiments of the invention are also directed to using recursive partitioning and decision trees to analyze genotype-phenotype associates and to identify putative strata with distinct patterns of interactions between genes and other non-genetic variables.

In embodiments of the invention, the methods comprise obtaining a sample containing genomic DNA from an individual for analysis. The sample can be, for example, a buccal swab, a saliva sample, a blood sample, or other suitable samples containing genomic DNA, as disclosed herein, and the like. The genomic DNA can be analyzed using any common technique available to the skilled person, such as, for example, high-throughput or low density array technologies, and the like. Results from such genotyping can subsequently be analyzed for the presence of certain variants known to be susceptibility variants for a particular condition, such as the genetic variants disclosed herein. Calculating risk conferred by a particular genotype for the individual can be based on comparing the genotype of the individual to previously determined risk (expressed as a relative risk (RR) or an odds ratio (OR), for example) for the genotype, for example for a heterozygous carrier of an at-risk variant for a particular condition or trait (such as for adverse effects from administered anesthesia). The calculated risk for the individual can be the relative risk for a person, or for a specific genotype of a person, compared to the average population with matched gender and ethnicity. The average population risk can be expressed as a weighted average of the risks of different genotypes, using results from a reference population, and the appropriate calculations to calculate the risk of a genotype group relative to the population can then be performed. Alternatively, the risk for an individual is based on a comparison of particular genotypes, for example heterozygous carriers of an at-risk allele of a marker compared with non-carriers of the at-risk allele. In some embodiments, using the population average can be more convenient, since it provides a measure which is easy to interpret for the user, i.e. a measure that gives the risk for the individual, based on his/her genotype, compared with the average in the population.

Overall risk for multiple risk variants can be performed using standard methodology. For example, assuming a multiplicative model, i.e. assuming that the risk of individual risk variants multiply to establish the overall effect, allows for a straight-forward calculation of the overall risk for multiple markers.

The detection of the particular genetic marker alleles that make up particular haplotypes in the sample can be performed by a variety of methods as described herein and/or known in the art. For example, genetic markers can be detected at the nucleic acid level (e.g. by direct nucleotide sequencing or by other genotyping means known to the skilled in the art) or at the amino acid level if the genetic marker affects the coding sequence of a protein (e.g. by protein sequencing or by immunoassays using antibodies that recognize such a protein). The marker alleles or haplotypes disclosed herein correspond to fragments of genomic segments (e.g. genes) associated with development of adverse effects of administered anesthetic, analgesia, and/or opioid. Such fragments encompass the DNA sequence of the polymorphic marker or haplotype in question but can also include DNA segments in strong LD (linkage disequilibrium) with the marker or haplotype. In some embodiments, such segments comprise segments in LD with the marker or haplotype as determined by a value of r² greater than 0.2 and/or |D′|>0.8).

In some embodiments, determination of susceptibility of developing adverse effects of administered anesthesia can be carried out using hybridization methods. (See Current Protocols in Molecular Biology, Ausubel, F. et al., eds., John Wiley & Sons, including all supplements). The presence of a specific marker allele can be indicated by sequence-specific hybridization of a nucleic acid probe specific for the particular allele. The presence of more than one specific marker allele or a specific haplotype can be indicated by using several sequence-specific nucleic acid probes, each being specific for a particular allele. A sequence-specific probe can be directed to hybridize to genomic DNA, RNA, or cDNA. A “nucleic acid probe” can be a DNA probe or an RNA probe that hybridizes to a complementary sequence. One of skill in the art would know how to design such a probe such that sequence specific hybridization will occur only if a particular allele is present in a genomic sequence from a test sample. The invention can also be reduced to practice using any convenient genotyping method, including commercially available technologies and methods for genotyping particular polymorphic markers.

A hybridization sample can be formed by contacting the test sample containing an anesthesia adverse effect-associated nucleic acid, such as a genomic DNA sample, with at least one nucleic acid probe. A non-limiting example of a probe for detecting mRNA or genomic DNA is a labeled nucleic acid probe that is capable of hybridizing to mRNA or genomic DNA sequences described herein. The nucleic acid probe can be, for example, a full-length nucleic acid molecule or a portion thereof, such as an oligonucleotide of at least 15, 30, 50, 100, 250, or 500 nucleotides in length that is sufficient to specifically hybridize under stringent conditions to appropriate mRNA or genomic DNA. In some embodiments, the nucleic acid probe can comprise at least one allele of at least one of the polymorphic markers selected from the group of: ABCB11 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-211C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs 8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, V158M rs4680, the CYP2D6 *2, *2A, *3, *4, *5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19,*20, *35, *40, *41, *42, and *44 alleles, the PGRN-Seq genetic markers (Table 16), and markers in linkage disequilibrium therewith, or the probe can be the complementary sequence of such a sequence. Hybridization can be performed by methods well-known to the person skilled in the art (see, e.g., Current Protocols in Molecular Biology, Ausubel, F. et al., eds., John Wiley & Sons, including all supplements). In some embodiments, hybridization refers to specific hybridization, i.e. hybridization with no mismatches (exact hybridization). In some embodiments, the hybridization conditions for specific hybridization are high stringency.

Specific hybridization, if present, is detected using standard methods. If specific hybridization occurs between the nucleic acid probe and the nucleic acid in the test sample, then the sample contains the allele that is complementary to the nucleotide that is present in the nucleic acid probe. The process can be repeated for any markers as disclosed herein, or markers that make up a haplotype as disclosed herein, or multiple probes can be used concurrently to detect more than one marker allele at a time. A single probe can also be designed in which the probe contains more than one marker allele of a particular haplotype (e.g. a probe containing alleles complementary to 2, 3, 4, 5, or all of the markers that make up a particular haplotype). Detection of the particular markers of the haplotype in the sample is indicative that the source of the sample has the particular haplotype (e.g. a haplotype) and therefore is susceptible or at risk of suffering adverse effects from administered anesthesia.

The detection probe can be of any suitable size, although preferably the probe is relatively short. In some embodiments, the probe is from 5-100 nucleotides in length. In some embodiments, the probe is from 10-50 nucleotides in length. In some embodiments, the probe is from 12-30 nucleotides in length. Other lengths of the probe are also contemplated and within the scope of the skill of the average person skilled in the art.

In some embodiments, the DNA template containing the SNP polymorphism is amplified by polymerase chain reaction (PCR) prior to detection. In such embodiments, the amplified DNA serves as the template for a detection probe and an enhancer probe.

In some embodiments, the detection probe, the enhancer probe, and/or the primers used for amplification of the template by PCR include the use of modified bases, including modified A and modified G. The use of modified bases can be useful for adjusting the melting temperature of the nucleotide molecule (probe and/or primer) to the template DNA, for example, for increasing the melting temperature in regions containing a low percentage of G or C bases, in which modified A with the capability of forming three hydrogen bonds to its complementary T can be used, or for decreasing the melting temperature in regions containing a high percentage of G or C bases, for example, by using modified G bases that form only two hydrogen bonds to their complementary C base in a double-stranded DNA molecule. In some embodiments, modified bases are used in the design of the detection nucleotide probe. Any modified base known to the skilled person can be selected in these methods, and the selection of suitable bases is well within the scope of the skilled person based on the teachings herein and known bases available from commercial sources as known to the skilled person.

In some embodiments, a peptide nucleic acid (PNA) probe can be used in addition to, or instead of, a nucleic acid probe in the hybridization methods described herein. A PNA is a DNA mimic having a peptide-like, inorganic backbone, such as N-(2-aminoethyl)glycine units, with an organic base (A, G, C, T, or U) attached to the glycine nitrogen via a methylene carbonyl linker (see, for example, Nielsen, P., et al., Bioconjug. Chem. 5:3-7 (1994)). The PNA probe can be designed to specifically hybridize to a molecule in a sample suspected of containing one or more of the marker alleles or haplotypes that are correlated with adverse effects of administered anesthesia. Hybridization of the PNA probe is thus diagnostic for susceptibility to such effects.

Embodiments of the invention are also directed to detecting SNPs within a set of genes, methods of detection include but are not limited to, for example, use of SNP microarrays, gene chips, dynamic allele-specific hybridization, molecular beacons, restriction fragment length polymorphism (RFLP)-based methods, PCR-based methods, flap endonuclease-based methods, primer extension, 5′-nuclease-based methods, oligonucleotide ligase assays, single strand conformation polymorphism methods, temperature gradient gel electrophoresis, denaturing high performance liquid chromatography, high-resolution melting analysis, DNA mismatch-binding methods, capillary electrophoresis, and next-generation sequencing methods, and the like.

In embodiments of the invention, a test sample containing genomic DNA obtained from the subject is collected, and PCR is used to amplify a fragment comprising one or more markers or haplotypes of the present invention. As disclosed herein, identification of a particular marker allele or haplotype can be accomplished using a variety of methods (e.g. sequence analysis, analysis by restriction digestion, specific hybridization, single-stranded conformation polymorphism assays (SSCP), electrophoretic analysis, and the like). In some embodiments, diagnosis is accomplished by expression analysis, for example by using quantitative PCR (kinetic thermal cycling). This technique can, for example, utilize commercially available technologies, such as TaqMan® (Applied Biosystems, Foster City, Calif.). The technique can assess the presence of an alteration in the expression or composition of a polypeptide or splicing variant(s). Further, the expression of the variant(s) can be quantified as physically or functionally different.

In some embodiments, the DNA template can be amplified by means of whole genome amplification (WGA) methods prior to assessment for the presence of specific polymorphic markers as described herein. Standard methods well-known to the skilled person for performing WGA can be utilized and are within scope of the invention.

In some embodiments, analysis by restriction digestion can be used to detect a particular allele if the allele results in the creation or elimination of a restriction site relative to a reference sequence. RFLP analysis can be conducted, e.g. as described in Current Protocols in Molecular Biology, supra. The digestion pattern of the relevant DNA fragment indicates the presence or absence of the particular allele in the sample.

Sequence analysis can also be used to detect specific alleles or haplotypes. Accordingly, in some embodiments, determination of the presence or absence of a particular marker allele or haplotype comprises sequence analysis of a test sample of DNA or RNA obtained from a subject or individual. PCR or other appropriate methods can be used to amplify a portion of a nucleic acid that contains a polymorphic marker or haplotype, and the presence of specific alleles can then be detected directly by sequencing the polymorphic site (or multiple polymorphic sites in a haplotype) of the genomic DNA in the sample.

In some embodiments, arrays of oligonucleotide probes that are complementary to target nucleic acid sequence segments from a subject can be used to identify particular alleles at polymorphic sites. For example, an oligonucleotide array can be used. Oligonucleotide arrays typically comprise a plurality of different oligonucleotide probes that are coupled to a surface of a substrate in different known locations. These arrays can generally be produced using mechanical synthesis methods or light-directed synthesis methods that incorporate a combination of photolithographic methods and solid phase oligonucleotide synthesis methods, or by other methods known to the person skilled in the art (see, e.g., Bier, et al. Adv Biochem Eng Biotechnol 109:433-53 (2008); Hoheisel, J. D., Nat Rev Genet. 7:200-10 (2006); Fan, et al. Methods Enzymol 410:57-73 (2006); Raqoussis, J. and Elvidge, G., Expert Rev Mol Diagn 6:145-52 (2006); and Mockler, et al. Genomics 85:1-15 (2005), each of which is incorporated herein by reference in its entirety). Many additional descriptions of the preparation and use of oligonucleotide arrays for detection of polymorphisms can be found, for example, in U.S. Pat. No. 6,858,394, U.S. Pat. No. 6,429,027, U.S. Pat. No. 5,445,934, U.S. Pat. No. 5,700,637, U.S. Pat. No. 5,744,305, U.S. Pat. No. 5,945,334, U.S. Pat. No. 6,054,270, U.S. Pat. No. 6,300,063, U.S. Pat. No. 6,733,977, U.S. Pat. No. 7,364,858, EP 619 321, and EP 373 203, each of which is incorporated herein by reference in its entirety.

Other methods of nucleic acid analysis that are available to those skilled in the art can be used to detect a particular allele at a polymorphic site. Representative methods include, for example, direct manual sequencing (Church and Gilbert, Proc. Natl. Acad. Sci. USA, 81: 1991-1995 (1988); Sanger, et al., Proc. Natl. Acad. Sci. USA, 74:5463-5467 (1977); Beavis, et al., U.S. Pat. No. 5,288,644); automated fluorescent sequencing; single-stranded conformation polymorphism assays (SSCPs); clamped denaturing gel electrophoresis (CDGE); denaturing gradient gel electrophoresis (DGGE) (Sheffield, et al., Proc. Natl. Acad. Sci. USA, 86:232-236 (1989)), mobility shift analysis (Orita, et al., Proc. Natl. AcadL Sci. USA, 86:2766-2770 (1989)), restriction enzyme analysis (Flavell, et al., Cell, 15:25-41 (1978); Geever, et al., Proc. Natl. AcadL Sci. USA, 78:5081-5085 (1981)); heteroduplex analysis; chemical mismatch cleavage (CMC) (Cotton, et al., Proc. Natl. Acad. Sci. USA, 85:4397-4401 (1985)); RNase protection assays (Myers, et al., Science, 230:1242-1246 (1985); use of polypeptides that recognize nucleotide mismatches, such as E. coli mutS protein; and allele-specific PCR, and the like. One of skill in the art can recognize other techniques that can be used for these purposes.

In embodiments of the invention, a determination of a susceptibility or risk of developing adverse effects to administered anesthesia, analgesic, and/or opioid can be made by examining expression and/or composition of a polypeptide encoded by a nucleic acid associated with adverse effects to administered anesthesia in those instances where the genetic marker(s) or haplotype(s) as disclosed herein result in a change in the composition or expression of the polypeptide. Thus, determination of a susceptibility to developing adverse effects to administered anesthesia, analgesic, and/or opioid can be made by examining expression and/or composition of one of these polypeptides or another polypeptide encoded by a nucleic acid associated with development of adverse effects to administered anesthesia, in those instances where the genetic marker or haplotype of the present invention results in a change in the composition or expression of the polypeptide.

A variety of methods can be used for detecting protein expression levels, including, for example, enzyme linked immunosorbent assays (ELISA), Western blots, immunoprecipitations, and immunofluorescence, and the like. A test sample from a subject can be assessed for the presence of an alteration in the expression and/or an alteration in composition of the polypeptide encoded by a particular nucleic acid. An alteration in the expression of a polypeptide encoded by the nucleic acid can be, for example, an alteration in the quantitative polypeptide expression (i.e. the amount of polypeptide produced). An alteration in the composition of a polypeptide encoded by the nucleic acid is an alteration in the qualitative polypeptide expression (e.g. expression of a mutant polypeptide or of a different splicing variant). In some embodiments, determination of a susceptibility to developing adverse effects to administered anesthesia is carried out by detecting a particular splicing variant encoded by a nucleic acid associated with development of such adverse effects or a particular pattern of splicing variants.

Both such alterations (quantitative and qualitative) can also be present. An “alteration” in the polypeptide expression or composition, as used herein, refers to an alteration in expression or composition in a test sample, as compared to the expression or composition of the polypeptide in a control sample. A control sample is a sample that corresponds to the test sample (e.g. is from the same type of cells) and is from a subject who is not affected by and/or who does not have a susceptibility to development of adverse effects to administered anesthesia. In some embodiments, the control sample is from a subject that does not possess a marker allele or haplotype associated with development of adverse effects to administered anesthesia, as disclosed herein. Similarly, the presence of one or more different splicing variants in the test sample or the presence of significantly different amounts of different splicing variants in the test sample, as compared with the control sample, can be indicative of a susceptibility to developing adverse effects to administered anesthesia. An alteration in the expression or composition of the polypeptide in the test sample, as compared with the control sample, can be indicative of a specific allele in the instance where the allele alters a splice site relative to the reference in the control sample. Various means of examining expression or composition of a polypeptide encoded by a nucleic acid are known to the person skilled in the art and can be used, including spectroscopy, colorimetry, electrophoresis, isoelectric focusing, and immunoassays (e.g. David et al., U.S. Pat. No. 4,376,110), such as immunoblotting (see, e.g., Current Protocols in Molecular Biology, particularly chapter 10, supra).

For example, in some embodiments, an antibody (e.g. an antibody with a detectable label) that is capable of binding to a polypeptide encoded by a nucleic acid associated with development of adverse effects to administered anesthesia can be used. Antibodies can be polyclonal or monoclonal. An intact antibody or a fragment thereof (e.g. Fv, Fab, Fab′, F(ab′)₂) can be used. The term “labeled,” with regard to the probe or antibody, encompasses direct labeling of the probe or antibody by coupling (i.e. physically linking) a detectable substance to the probe or antibody, as well as indirect labeling of the probe or antibody by reactivity with another reagent that is directly labeled. Examples of indirect labeling include detection of a primary antibody using a labeled secondary antibody (e.g. a fluorescently-labeled secondary antibody) and end-labeling of a DNA probe with biotin such that it can be detected with fluorescently-labeled streptavidin.

In some embodiments, the level or amount of a polypeptide in a test sample is compared with the level or amount of the polypeptide in a control sample. A level or amount of the polypeptide in the test sample that is higher or lower than the level or amount of the polypeptide in the control sample, such that the difference is statistically significant, is indicative of an alteration in the expression of the polypeptide encoded by the nucleic acid and is diagnostic for a particular allele or haplotype responsible for causing the difference in expression. Alternatively, the composition of the polypeptide in a test sample is compared with the composition of the polypeptide in a control sample. In some embodiments, both the level or amount and the composition of the polypeptide can be assessed in the test sample and in the control sample.

In some embodiments, determination of a susceptibility to developing adverse effects to administered anesthesia is made by detecting at least one marker or haplotype as disclosed herein, in combination with an additional protein-based, RNA-based, or DNA-based assay. In some embodiments, determination of a susceptibility to inadequate pain relief or sensitivity is made by detecting at least one marker or haplotype as disclosed herein, in combination with an additional protein-based, RNA-based, or DNA-based assay. In some embodiments, determination of a susceptibility to OSA is made by detecting at least one marker or haplotype as disclosed herein, in combination with an additional protein-based, RNA-based, or DNA-based assay. In some embodiments, determination of a susceptibility to developing adverse effects to administered anesthesia, analgesic, and/or opioid is made by determining incidence of OSA made by detecting at least one marker or haplotype as disclosed herein, in combination with an additional protein-based, RNA-based, or DNA-based assay.

In some embodiments, determination of a susceptibility to developing adverse effects to an administered opioid is made by detecting at least one marker or haplotype as disclosed herein, in combination with an additional protein-based, RNA-based, or DNA-based assay. In some embodiments, the administered opioid is morphine. In some embodiments, the administered opioid is codeine. In some embodiments, characterization of the OCT1 genotype is used to determine susceptibility to developing adverse effects to an administered opioid. In some embodiments, the UGT2B7 rs7438135 allele is determined. In some embodiments, the rs622342 allele is characterized.

Kits

Kits useful in the methods of the invention comprise components useful in any of the methods disclosed herein, including, for example, primers for nucleic acid amplification, hybridization probes, restriction enzymes (e.g. for RFLP analysis), allele-specific oligonucleotides, antibodies that bind to an altered polypeptide encoded by a nucleic acid as disclosed herein (e.g. a genomic segment comprising at least one polymorphic marker and/or haplotype of the present invention) or to a non-altered (native) polypeptide encoded by a nucleic acid of the invention as described herein, means for amplification of a nucleic acid associated with the development of adverse effects to administered anesthesia, means for analyzing the nucleic acid sequence of a nucleic acid associated with development of adverse effects to administered anesthesia, means for analyzing the amino acid sequence of a polypeptide encoded by a nucleic acid associated with development of adverse effects to administered anesthesia, and the like. The kits can, for example, include necessary buffers, nucleic acid primers for amplifying nucleic acids of the invention (e.g. a nucleic acid segment comprising one or more of the polymorphic markers as described herein), and reagents for allele-specific detection of the fragments amplified using such primers and necessary enzymes (e.g. DNA polymerase), and the like. Additionally, kits can provide reagents for assays to be used in combination with the methods as disclosed herein, e.g. reagents for use with other diagnostic assays for determining susceptibility to development of adverse effects to administered anesthesia. In some embodiments, reagents for performing WGA are included in the reagent kit.

In some embodiments, a kit for assaying a sample from a subject to detect a risk or susceptibility to developing adverse effects to administered anesthesia, analgesic, and/or opioid in a subject is provided, wherein the kit comprises reagents necessary for selectively detecting at least one allele of at least one polymorphism as disclosed herein in the genome of the individual. In some embodiments, the reagents comprise at least one contiguous oligonucleotide that hybridizes to a fragment of the genome of the individual comprising at least one polymorphism, as disclosed herein. In some embodiments, the reagents comprise at least one pair of oligonucleotides that hybridize to opposite strands of a genomic segment obtained from a subject, wherein each oligonucleotide primer pair is designed to selectively amplify a fragment of the genome of the individual that includes at least one polymorphism associated with developing adverse effects to administered anesthesia. In some embodiments, reagents comprise at least one pair of oligonucleotides that hybridize to opposite strands of a genomic segment obtained from a subject, wherein each oligonucleotide primer pair is designed to selectively amplify a fragment of the genome of the individual that includes at least one polymorphism associated with developing adverse effects to an administered opioid. In some embodiments, the polymorphism is selected from: ABCB1 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-211C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, V158M rs4680, the CYP2D6 *2, *2A, *3, *4, *5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19, *20, *35, *40, *41, *42, and *44 alleles, the PGRN-Seq genetic markers (Table 16), and markers in linkage disequilibrium therewith. In some embodiments, the fragment is at least 20 base pairs in size. Such oligonucleotides or nucleic acids (e.g. oligonucleotide primers) can be designed using portions of the nucleic acid sequence flanking polymorphisms (e.g. SNPs or microsatellites) that are associated with a risk of developing adverse effects to administered anesthesia. In some embodiments, the kit comprises one or more labeled nucleic acids capable of allele-specific detection of one or more specific polymorphic markers or haplotypes, as well as reagents for detection of the label. Suitable labels include, e.g. a radioisotope, a fluorescent label, an enzyme label, an enzyme co-factor label, a magnetic label, a spin label, an epitope label, and the like.

In some embodiments, the polymorphic marker or haplotype to be detected by the reagents of the kit comprises one or more markers, two or more markers, three or more markers, four or more markers, or five or more markers selected from the group consisting of: ABCB1 rs1045642 (ABCB1 C3435T), ABCB1 G2677T/A, ABCC3-211C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, the CYP2D6 *2, *2A, *3, *4, *5, *6, *7, *8, *9, *10, *11, *14, *15, *17, *18, *19, *20, *35, *40,*41, *42, and *44 alleles, the PGRN-Seq genetic markers (Table 16), and markers in linkage disequilibrium therewith. In some embodiments, the marker or haplotype to be detected comprises one or more markers, two or more markers, three or more markers, four or more markers, or five or more markers selected from the group consisting of: the markers ABCB11 rs1045642 (ABCB11 C3435T), ABCB1 G2677T/A, ABCC3-211C>T rs4793665, ADRB2 rs1042717, ADRB2 rs1042714, ADRB2 rs1042713, ANKKI rs1800497, DRD2 rs6279, FAAH rs4141964, FAAH rs2295632, FAAH rs3766246, FAAH rs324420, FAAH rs932816, FAAH rs324419, COMT rs6269, COMT rs4818, COMT rs4680 (COMT Val158Met), COMT rs4633, GCHI rs 8007267, GCHI rs752688, GCHI rs4411417, OPRM1 A118G, CYP2D6, MDR1 rs1045642, MDR1 rs9282564, MDR1 rs1128503, MDR1 rs2032582, MDR1 rs2229109, OCT1 rs12208357, OCT1 rs34130495, OCT1 rs72552763, OCT1 rs34059508, TRPA1 rs1443952, TRPA1 rs13279503, TRPA1 rs13255063, TRPA1 rs1947913, UGT2B7 rs7439366, UGT2B7 rs7668258, V158M rs4680, the CYP2D6 *2, *2A, *3, *4, *5, *6, *7, *8, *9, * 10,*11, *14, *15, *17, *18, *19, *20, *35, *40, *41, *42, and *44 alleles, the PGRN-Seq genetic markers (Table 16), and markers in linkage disequilibrium therewith.

In certain embodiments, the kit further comprises a set of instructions for using the reagents comprising the kit.

Anesthetics and Analgesics

Embodiments of the invention are also directed to assessing patient response to anesthesia, anesthetics not limited to local anesthetics, general anesthetics, inhaled agents, intravenous agents, and muscle relaxants.

“Administered” anesthesia includes any anesthetic that is introduced to a subject, including, but not limited to, anesthesia that is administered orally, by dermal contact, subcutaneously, intravenously, by epidural means, spinally, and by inhalation.

Local anesthetics include, but are not limited to, benzocaine, chloroprocaine, cocaine, cyclomethycaine, dimethocaine/loracaine, piperocaine, propoxycaine, procaine/novocaine, proparacaine, tetracaine/amethocaine, articaine, bupivacaine, cinchocaine/dibucaine, etidocaine, levobupivacaine, lidocaine/lignocaine, mepivacaine, prilocaine, ropivacaine, trimecaine, saxitoxin, tetrodotoxin, menthol, eugenol, and combinations and derivations thereof.

General anesthetics include, but are not limited to, desflurane, enflurane, halothane, isoflurane, methoxyflurane, nitrous oxide, sevoflurane, xenon, amobarbital, methohexital, thiamylal, thiopental, diazepam, lorazepam, midazolam, etomidate, ketamine, propofol, alfentanil, fentanyl, remifentanil, sufentanil, buprenorphine, butorphanol, diamorphine, hydromorphone, levorphanol, meperidine, methadone, morphine, nalbuphine, oxycodone, oxymorphone, pentazocine, succinylcholine, decamethonium, mivacurium, rapacuronium, atracurium, cisatracurium, rocuronium, vecuronium, alcuronium, doxacurium, gallamine, metocurine, pancuronium, pipecuronium, tubocurarine, and combinations and derivations thereof.

Analgesics include, but are not limited to, opioids, such as morphine, fentanyl, hydromorphone, codeine, oxycodone, hydrocodone, tramadol, ondansetron, dexamethasone, methadone, alfentanil, remifentanil, and derivations thereof, and non-opioid analgesics, including acetaminophen, NSAIDs, dexmedetomidine, clonidine, and combinations and derivations thereof.

Adverse Anesthetic and Analgesic Effects

Adverse anesthetic effects can include, but are not limited to, inadequate pain relief, increased pain sensitivity, increased anesthesia and analgesic requirement, PONV, respiratory depression, excessive sedation, pruritus, death, and the like.

Genetic Indication for Codeine-Related Adverse Effects and Analgesia

As described herein, high CYP2D6 activity score was associated with increased adverse drug reactions (ADRs) (p=0.004) during post-operative days (POD) 0-2. Sedation after codeine was more common in girls (p=0.05). High pain intensity (p=0.003) and an interaction between CYP2D6 activity and time of the day after surgery (p=0.003) contributed to sedation risk. Pain reduction following a dose of codeine was associated with age (p=0.0002) and time (p=0.001) but not CYP2D6.

The results described herein, together with previous reports of codeine-related deaths and serious ADRs, demonstrate CYP2D6 association with codeine-related ADRs and sedation and indicate the potential safety hazards of using unpredictable and potentially life-threatening opioids that are also metabolized through the CYP2D6 pathway, such as codeine, hydrocodone, oxycodone, and tramadol. CYP2D6 testing prior to opioid administration is necessary to avoid codeine-, hydrocodone-, oxycodone-, or tramadol-related death and severe ADR risks, especially in young children and infants of breastfeeding mothers who have increased metabolic capacity to convert these opioids to their more potent metabolites. These should not be considered safe alternatives without CYP2D6 testing. An alternative approach to CYP2D6 testing is prescribing relatively safer non-opioid analgesics such as acetaminophen, non-steroidal anti-inflammatory drugs or appropriate and as needed doses of oral morphine in high-risk pediatric population.

Accordingly, embodiments of the invention are also directed to assessing a patient's risk of adverse response to opioids. In some embodiments, the patient's risk of adverse response to an opioid is assessed by CYP2D6 testing. In some embodiments, the opioid is codeine. In some embodiments, the opioid is, for example, hydrocodone, oxycodone, tramadol, or the like, or other commonly used medications used to prevent postoperative nausea and vomiting, such as, for example, ondansetron and dexamethasone, and the like.

In some embodiments, the patient at risk for an adverse response to an opioid is a child. In some embodiments, the patient at risk for an adverse response to an opioid is an infant.

Embodiments of the invention are also directed to adapting a treatment strategy for a patient requiring pain relief. In some embodiments, the patient's risk of adverse response to an opioid is assessed by CYP2D6 testing, and results of the risk assessment are used to make a decision regarding administration of an opioid. In some embodiments, an opioid is administered to a patient with low risk of adverse response to an opioid according to CYP2D6 testing. In some embodiments, an alternative therapeutic is administered to a patient with high risk of adverse response to an opioid according to CYP2D6 testing.

Having described the invention in detail, it will be apparent that modifications, variations, and equivalent embodiments are possible without departing from the scope of the invention defined in the appended claims. All examples in the present disclosure are provided as non-limiting examples.

EXAMPLES

The following non-limiting examples are provided to further illustrate embodiments of the invention disclosed herein. It will be appreciated by those of skill in the art that the techniques disclosed in the examples that follow represent approaches that have been found to function well in the practice of the invention and thus can be considered to constitute examples of modes for its practice. However, those of skill in the art will, in light of the present disclosure, appreciate that many changes can be made in the specific embodiments that are disclosed and still obtain a like or similar result without departing from the spirit and scope of the invention.

Example 1 Genotypes Associated with Opioid-Related Adverse Effects

Individuals can have highly variable responses to opioids (FIG. 2). A perioperative morphine pharmacogenetic study of 150 children undergoing tonsillectomy revealed associations between genetic/non-genetic factors and postoperative opioid-related adverse effects and inadequate analgesia. The mean age of study participants was 9.25±2.65 years, the ratio of boys to girls was 47:53, the mean weight was 39.3±15.1 kg, the race ratio of white to non-white children was 81%:19%, and the ratio of children with obstructive sleep apnea (OSA) to those with no OSA was 46:54.

After adjusting for non-genetic risk factors such as OSA, the TI genotype of the ABCB1 SNP (rs1045642) was associated with higher risks of respiratory depression (OR=2.2, p=0.035) and PONV (OR=3.45, p=0.029) from morphine compared to the CT genotype (See Table 1 and FIG. 3). The CC genotype of ABCB1 was protective against respiratory depression (OR=0.46) and PONV (OR=0.29). The AA genotype of fatty acid amide hydrolase (FAAH) (rs4141964) is associated higher incidence of PONV (OR: 3.55, p=0.0164) compared to the GA genotype, and the GG genotype of FAAH is protective against PONV (See Table 2 and FIG. 4). There was an over ten-fold difference in the occurrence of opioid-induced PONV between high risk genotypes (IT of ABCB1 (rs1045642); AA of FAAH (rs4141964)) and low risk genotypes (CC of ABCB1 and GG of FAAH) (FIG. 4).

TABLE 1 PACU Incidence of Prolonged PACU Genotype Opioid Side Effects (N, %) Stay (N, proportion) ABCB1 Respiratory Excessive Respiratory Excessive rs1045642 Depression Sedation PONV Depression Sedation PONV CC  9 (21%) 7 (18%) 7 (17%) 2 (5%) 12 (28%) 0 (0%) CT 17 (26%) 13 (20%)  7 (11%) 9 (14%) 22 (33%) 4 (6%) TT 12 (29%) 5 (14%) 8 (20%) 7 (17%) 13 (32%) 5 (12%)

TABLE 2 PACU Incidence Prolonged PACU Genotype of Opioid Side Effects (N, %) Stay (N, proportion) FAAH Excessive Respiratory Excessive rs4141964 Respiratory Depression Sedation PONV Depression Sedation PONV GG 9 (16%) 6 (11%) 4 (7%) 3 (5%) 12 (21%) 0 (0%) GA 20 (31%)  16 (27%)  12 (19%) 9 (14%) 26 (41%) 5 (8%) AA 9 (30%) 3 (11%)  6 (20%) 6 (20%)  9 (30%) 4 (13%)

Example 2 Genotypic Association with Opioid-Induced Respiratory Depression

In the same perioperative morphine pharmacogenetic study (Example 1), it was observed that after exposure to morphine, the resting minute ventilation (MV) decreased by 47.5% in children with the ABCB1 T genotype, compared to only 18.4% in those with the CC and CT genotypes (FIG. 5). The trend persisted after a 5% carbon dioxide challenge, as resting MV after morphine decreased 29% in children with the ABCB1 T genotype, compared to 10% in those with the CC and CT genotypes (p<0.05) (FIG. 5). This further illustrates that the presence of the TT genotype of ABCB1 (rs1045642) is associated with development of adverse effects from administered anesthesia.

Example 3 Gene-Gene Interactions and Increased Risk of Opioid-Related Respiratory Depression

In the same perioperative morphine pharmacogenetic study (Example 1), it was observed that, for children suffering from prolonged respiratory depression in the PACU, a significant interaction was detected between ABCB1 SNP (rs1045642) and FAAH SNP (rs4141964) (p=0.022). The probability of having the highest risk of respiratory depression (46.7%) is observed in the combination of TI and AA genotypes in ABCB1 (rs1045642) and FAAH (rs4141964), respectively (FIG. 6). The lowest probability of respiratory depression (1.5%) was observed in the combination of CC and AA genotypes in ABCB1 (rs1045642) and FAAH (rs4141964), respectively (p=0.022) (FIG. 6).

Example 4 Economic Burden of a Prolonged PACU Stay

The AA genotype of FAAH (rs4141964) was associated with longer PACU stays (143.8±14 minutes, or approximately USD $570), compared to GA (125.6±12.7 minutes, or approximately USD $480) and GG (112.5±10.8 minutes, or approximately USD $410) (p=0.0178). On average, children with AA genotypes stayed 24 minutes longer in the PACU than those with GG genotypes. This illustrates the need for creating and/or developing personalized anesthesia protocols in order to reduce development of adverse effects to administered anesthesia and consequently, hospital costs associated with caring for patients who suffer from such effects.

Example 5 Prolonged PACU Stays

In the same perioperative morphine pharmacogenetic study (Example 1), four SNPs of COMT had strong associations with pain scores, defined as maximum face, legs, activity, cry, and consolability (FLACC), morphine requirement, and analgesic interventions in the PACU (see Table 3).

TABLE 3 ^(a)Postoperative ^(b)Total ^(c)Postoperative Pain Morphine Use Opioid COMT Genotype (FLACC) (mg/kg) Intervention SNP Frequency mean □ SE p-value mean □ SE p-value OR (95% CI) p-value rs6269 AA(49) 2.1 □ 0.4 0.0205 0.23 □ 0.01 0.0124 0.30 (0.14-0.64) 0.0017 AG + GG(101) 3.2 □ 0.3 0.25 □ 0.01 rs4818 CC (58) 2.4 □ 0.4 0.1356 0.23 ± 0.01 0.1193 0.39 (0.19-0.79) 0.0097 CG + GG (92) 3.1 □ 0.3 0.25 ± 0.01 rs4680 AA(38) 2.2 □ 0.4 0.1139 0.23 □ 0.01 0.0568 0.32 (0.15-0.72) 0.0053 AG + GG(112) 3.1 □ 0.3 0.25 □ 0.01 rs4633 TT (37) 2.1 □ 0.4 0.0405 0.22 □ 0.01 0.0134 0.30 (0.13-0.67) 0.0032 TC + CC (113) 3.1 □ 0.3 0.25 □ 0.01 ^(a)adjusted for age, race and intra-operative morphine requirement ^(b)adjusted for race and intra-operative morphine requirement ^(c)adjusted for intra-operative morphine requirement

Following a standard intraoperative morphine dose, patients with certain genotypes of COMT SNPs (AG and GG of rs6269; CG and GG of rs4818; AG and GG of rs4680 and TC and CC of rs4633) had higher surgical pain and were 3 times more likely to require analgesic interventions compared to other respective genotypes in the PACU (Table 3).

Example 6 Additional Genotypes Associated with Surgical Pain, Analgesia, and Opioid-Related Adverse Effects

In the same study, but with a larger sample size, additional genotypes were subsequently identified for further stratification of African-American and Caucasian children. These genotypes can be associated with specific phenotypes, such as, for example, FLAAC scores, inadequate pain control with long PACU stays, need for opioid analgesic intervention in the PACU, need for post-operative analgesia in the PACU, opioid-related side effects leading to prolonged PACU stays, opioid-related side effects, respiratory depression, respiratory depression resulting in prolonged PACU stays, and post-operative vomiting resulting in prolonged PACU stays.

As shown in Table 4, a number of genotypes were strongly associated with outcomes. These include the T allele of UGT2B7 rs7439366 (A/T), which was strongly associated with maximum PACU FLACC scores for African-American children. In addition, the C allele of TRPA1 rs13279503 (C/G), the A allele of TRPA1 rs13255063 (A/T), and the G allele of MDR1 (ABCB1) rs9282564 (G/A) were strongly associated with inadequate pain control with long PACU stays for African-American children.

The G allele of COMT rs4818 (G/C) and the A allele of FAAH rs324419 (A/G) were strongly associated with a need for opioid analgesic intervention in the PACU for African-American children. The G allele of MDR1(ABCB1) rs9282564 (G/A) was also associated with opioid-related side effects leading to prolonged PACU stays for African-American children, as was the A allele of FAAH rs932816 (A/G). The A allele of MDR1 (ABCB1) rs2229109 (A/C) was strongly associated with opioid-related side effects for Caucasian children.

The T allele of GCHI rs8007268 (T/C) and the A allele of ADRB2 rs1042717 (A/G) were both strongly associated with respiratory depression for African-American children. TRPA1 rs1947913 (A/T), MDR1 (ABCB1) rs9282564 (G/A), and FAAH rs932816 (A/G) were all strongly associated with respiratory depression resulting in prolonged PACU stays. Of these, the A allele of TRPA rs1947913 was a risk factor for Caucasian children, the A allele of FAAH rs932816 was a risk factor for African-American children, and the G allele of MDR1(ABCB1) rs9282564 was a risk factor for both Caucasian and African-American children.

Four FAAH SNPs, namely rs2295632 (A/C), rs4141964 (A/G), rs3766246 (T/C), and rs324420 (A/C), were strongly associated with post-operative vomiting resulting in prolonged PACU stays. The A, A, T, and A alleles, respectively, were determined to be risk factors for Caucasian children.

TABLE 4 Minor/major Risk Effect Outcomes SNP Gene Cohort P-value allele allele size Effect measure Maximum PACU rs1042713 ADRB2 Caucasian 0.034 A/G A 1.200934 ratio of counts FLACC scores rs2032582 MDR1 (ABCB1) Caucasian 0.037 AT/G G 0.837696 ratio of counts rs8007267 GCH1 African- 0.035 T/C C 0.817258 ratio of counts American rs1042713 ADRB2 interaction 0.029 A/G rs1042714 ADRB2 combined 0.046 G/C C 0.865455 ratio of counts Maximum PACU rs7439366 UGT2B7 African- 0.011 A/T T 11.36002 odds ratio American FLACC scores rs1800497 ANKK1 Caucasian 0.046 T/C T 0.470246 odds ratio zero-inflated rs7439366 UGT2B7 interaction 0.032 A/T Maximum PACU rs2295632 FAAH Caucasian 0.034 A/C A 0.6723 mean difference numerical rating rs4141964 FAAH Caucasian 0.05 A/G A 0.5889 mean difference scale rs3766246 FAAH Caucasian 0.035 T/C T 0.6336 mean difference rs1042713 ADRB2 Caucasian 0.046 A/G A 0.6577 mean difference rs2032582 MDR1 (ABCB1) African- 0.006 AT/G G −3.5215 mean difference American rs4680 COMT African- 0.032 A/G G −2.147 mean difference American rs8007267 GCH1 African- 0.043 T/C T 1.6818 mean difference American rs2032582 MDR1 (ABCB1) interaction 0.004 AT/G rs4680 COMT interaction 0.018 A/G rs324420 FAAH combined 0.029 A/C A 0.6744 mean difference Inadequate pain rs1042713 ADRB2 Caucasian 0.022 A/G A 1.648227 odds ratio control with long rs13279503 TRPA1 African- 0.034 C/G C 3.43941 odds ratio PACU stay American rs13255063 TRPA1 African- 0.008 A/T A 14.13424 odds ratio American rs9282564 MDR1 (ABCB1) African- 0.029 G/A G 9.042177 odds ratio American rs13255063 TRPA1 interaction 0.023 A/T rs9282564 MDR1 (ABCB1) interaction 0.039 G/A Opioid analgesic rs1042713 ADRB2 Caucasian 0.032 A/G A 1.575701 odds ratio intervention need rs4818 COMT African- 0.003 G/C G 14.78333 odds ratio in PACU American rs6279 DRD2 African- 0.027 C/G G 0.120308 odds ratio American rs324419 FAAH African- 0.006 A/G A 43.0215 odds ratio American rs4818 COMT interaction 0.008 G/C rs1800497 ANKK1 combined 0.037 T/C T 1.614782 odds ratio Postop analgesia rs6269 COMT Caucasian 0.022 G/A G 1.280051 count ratio in PACU Postoperative rs4141964 FAAH Caucasian 0.015 A/G A 0.0146 mean difference morphine use rs3766246 FAAH Caucasian 0.008 T/C T 0.0158 mean difference rs6269 COMT Caucasian 0.03 G/A G 0.0135 mean difference rs4633 COMT Caucasian 0.033 T/C C −0.0135 mean difference rs4818 COMT African- 4E−04 G/C C −0.0477 mean difference American rs4818 COMT interaction <0.0001 G/C rs6269 COMT interaction 0.014 G/A rs4141964 FAAH interaction 0.04 A/G rs2229109 MDR1 (ABCB1) combined 0.05 A/C C −0.0354 mean difference Opioid side rs2295632 FAAH Caucasian 0.042 A/C A 1.690966 odds ratio effects leading to rs4141964 FAAH Caucasian 0.028 A/G A 1.739504 odds ratio prolonged PACU rs3766246 FAAH Caucasian 0.031 T/C A 1.723919 odds ratio stay rs6269 COMT African- 0.038 G/A A 0.182246 odds ratio American rs9282564 MDR1 (ABCB1) African- 0.049 G/A G 12.9061 odds ratio American rs932816 FAAH African- 0.039 A/G A 6.073287 odds ratio American rs932816 FAAH interaction 0.012 A/G rs9282564 MDR1 (ABCB1) combined 0.015 G/A G 2.429535 odds ratio Opioid-related rs2229109 MDR1(ABCB1) Caucasian 0.034 A/C A 5.401086 odds ratio side effects rs7668258 UGT2B7 Caucasian 0.022 T/C T 1.918992 odds ratio Respiratory rs2295632 FAAH Caucasian 0.034 A/C A 1.764205 odds ratio depression rs1128503 MDR1(ABCB1) Caucasian 0.033 T/C C 0.556549 odds ratio rs8007267 GCH1 African- 0.006 T/C T 6.621355 odds ratio American rs1042713 ADRB2 African- 0.044 A/G G 0.379538 odds ratio American rs1042717 ADRB2 African- 0.036 A/G A 2.823847 odds ratio American rs8007267 GCH1 interaction 0.002 T/C rs1042713 ADRB2 interaction 0.017 A/G rs1128503 MDR1 (ABCB1) interaction 0.029 T/C rs1042717 ADRB2 interaction 0.048 A/G rs752688 GCH1 combined 0.029 T/C C 0.493417 odds ratio rs4411417 GCH1 combined 0.029 C/T T 0.493417 odds ratio Respiratory rs13279503 TRPA1 Caucasian 0.046 C/G C 1.842273 odds ratio depression rs1947913 TRPA1 Caucasian 0.033 A/T A 2.053612 odds ratio resulting in rs9282564 MDR1 (ABCB1) Caucasian 0.002 G/A G 3.724379 odds ratio prolonged PACU rs9282564 MDR1 (ABCB1) African- 0.027 G/A G 18.02032 odds ratio stay American rs932816 FAAH African- 0.043 A/G A 5.743679 odds ratio American rs932816 FAAH interaction 0.016 A/G rs1947913 TRPA1 interaction 0.017 A/T rs324420 FAAH interaction 0.034 A/C Postoperative rs7668258 UGT2B7 Caucasian 0.041 T/C C 0.493072 odds ratio vomiting resulting rs2295632 FAAH Caucasian 0.004 A/C A 2.578738 odds ratio in prolonged rs4141964 FAAH Caucasian 0.02 A/G A 2.177767 odds ratio PACU stay rs3766246 FAAH Caucasian 0.021 T/C T 2.168856 odds ratio rs324420 FAAH Caucasian 0.039 A/C A 2.021824 Significance using p < 0.05 

1-20. (canceled)
 21. A method for treating a human patient for post-operative pain, the method comprising the steps of a) obtaining a biological sample from a patient in need of treatment for post-operative pain wherein the biological sample is obtained prior to administering an opioid analgesic to the patient, b) genotyping the biological sample from the patient to determine the patient's genotype at a single nucleotide polymorphism (SNP) of the ATP binding cassette C3 (ABCC3) gene identified by rs4793665 or an SNP of the mu opioid receptor (OPRM1) gene designated A118G, wherein the genotyping is performed using one or more of a PCR-based technique for nucleic acid amplification from the sample, gel electrophoresis, mass spectrometry, and microarray hybridization, c) identifying the presence of a T allele at ABCC3 rs4793665 in the biological sample or d) identifying the presence of a G allele at OPRM1 A118G in the biological sample, and e) administering an analgesic selected from fentanyl, hydromorphone, codeine, oxycodone, hydrocodone, tramadol, ondansetron, dexamethasone, methadone, alfentanil, remifentanil, acetaminophen, a non-steroidal anti-inflammatory drug (NSAID), dexmedetomidine, clonidine, and combinations and derivations thereof for post-operative pain to the patient having an allele identified in step c) or f) administering a non-opioid analgesic for post-operative pain to the patient having an allele identified in step d).
 22. A method for treating a human patient for post-operative pain, the method comprising the steps of a) obtaining a biological sample from a patient in need of treatment for post-operative pain wherein the biological sample is obtained prior to administering an opioid analgesic to the patient, b) genotyping the biological sample from the patient to determine the patient's genotype at a single nucleotide polymorphism (SNP) of the ATP binding cassette B1 (ABCB1) gene identified by rs1045642 or an SNP of the fatty acid amide hydrolase (FAAH) gene identified by rs2295632, wherein the genotyping is performed using one or more of a PCR-based technique for nucleic acid amplification from the sample, gel electrophoresis, mass spectrometry, and microarray hybridization, c) identifying the presence of a T allele at ABCB1 rs1045642 in the biological sample or d) identifying the presence of an A allele at FAAH rs2295632 in the biological sample, and e) administering hydromorphone for post-operative pain to the patient having an allele identified in step c) or f) administering a non-opioid analgesic for post-operative pain to the patient having an allele identified in step d).
 23. A method for treating a human patient for post-operative pain, the method comprising the steps of a) obtaining a biological sample from a patient in need of treatment for post-operative pain wherein the biological sample is obtained prior to administering an opioid analgesic to the patient, b) genotyping the biological sample from the patient to determine the patient's genotype at a single nucleotide polymorphism (SNP) of the organic cation transporter 1 (OCT1) gene identified by rs12208357, wherein the genotyping is performed using one or more of a PCR-based technique for nucleic acid amplification from the sample, gel electrophoresis, mass spectrometry, and microarray hybridization, c) identifying the presence of a T allele at OCT1 rs12208357 in the biological sample, and d) administering fentanyl, hydromorphone, codeine, oxycodone, hydrocodone, tramadol, ondansetron, dexamethasone, methadone, alfentanil, remifentanil, acetaminophen, a non-steroidal anti-inflammatory drug (NSAID), dexmedetomidine, clonidine, and combinations and derivations thereof for post-operative pain to the patient having an allele identified in step c).
 24. The method of claim 23, further comprising genotyping the biological sample from the patient to determine the patient's genotype at one or more SNPs of the OCT1 gene identified by rs34130495, rs72552763, and rs34059508.
 25. The method of claim 21, wherein the patient is a pediatric patient.
 26. The method of claim 22, wherein the patient is a pediatric patient.
 27. The method of claim 23, wherein the patient is a pediatric patient.
 28. The method of claim 21, wherein the non-opioid analgesic is selected from the group consisting of acetaminophen, a non-steroidal anti-inflammatory drug (NSAID), dexmedetomidine, clonidine, and combinations and derivations thereof.
 29. The method of claim 22, wherein the non-opioid analgesic is selected from the group consisting of acetaminophen, a non-steroidal anti-inflammatory drug (NSAID), dexmedetomidine, clonidine, and combinations and derivations thereof.
 30. The method of claim 21, further comprising receiving at least one non-genetic risk factor selected from the patient's race, sex, and obstructive sleep apnea disease status.
 31. The method of claim 22, further comprising receiving at least one non-genetic risk factor selected from the patient's race, sex, and obstructive sleep apnea disease status.
 32. The method of claim 23, further comprising receiving at least one non-genetic risk factor selected from the patient's race, sex, and obstructive sleep apnea disease status.
 33. The method of claim 30, wherein the non-genetic risk factor is race and the race is selected from Caucasian and African American.
 34. The method of claim 33, wherein the race is Caucasian.
 35. The method of claim 31, wherein the non-genetic risk factor is race and the race is selected from Caucasian and African American.
 36. The method of claim 35, wherein the race is Caucasian.
 37. The method of claim 32, wherein the non-genetic risk factor is race and the race is selected from Caucasian and African American.
 38. The method of claim 37, wherein the race is Caucasian.
 39. The method of claim 30, wherein the non-genetic risk factor is obstructive sleep apnea disease status.
 40. The method of claim 31, wherein the non-genetic risk factor is obstructive sleep apnea disease status. 